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ABSTRACT 

This dissertation is an integrated method of designing sampled- 
data feedback control systems to required time response, This method 
is distinguished from the classical frequency design methods in two 
ways - first, the closed-loop rather than the open-loop function serves 
as the starting point of the method; and, second, the procedure invol- 
ves nearly direct synthesis rather than trialand error, To implement 
the design procedures, a set of correlation theorems are derived which 
permit interpretation of the closed-loop system time responses in 
terms of the pole-zero configuration of the closed-loop system function, 
Further, a graphical method is presented which permits rapid solu- 
tion of the open-loop poles when the closed-loop poles are known, The 
techniques and results of an analog computer study are presented which 
serve as an experimental verification of the methods developed in this 


dissertation, 


il 





IE 


1 


IV 


TABLE OF CONTENTS 


INTRODUCTION 

1,1 General Background 

1.2 Discussion of the Dissertation 
INTRODUCTION TO SAMPLED-DATA THEORY 
2.1 Operation of the Sampler 

2.2 Elements of z-Transform Theory 

2.3 z-Transfer Function 


2.4 Final and Initial Value Theorems, Error Coefficients 
and Mapping 


2.5 Response between the Sampling Instants 


2.6 Analysis and Design of Sampled-Data Control System 
by Conventional Techniques 


SYNTHESIS THROUGH CLOSED-LOOP TRANSFER 
FUNCTION 


3.1 Background Discussion 


3,2 Synthesis through the Closed-Loop Poles-Zeros 
Configuration 


3,3 Time Response and Basic Pole-Zero Configuration 
SYNTHESIS TECHNIQUES 

4.1 Introduction 

4,2 Correlation Theorems 

4,3 The Required Number of Poles and Zeros 

4,4 Open- Loop Root Locus 

4.5 Closed-Loop Root Locus 

4,6 Additional Design Aids 

APPLICATION OF SYNTHESIS PROCEDURE 


5.1 Introduction 


Page 


10 
23 


27 


34 


36 


40 


48 


50 


54 


58 
59 
77 
81 
90 
91 


98 





Page 





5,2 Complete Example of Synthesis Procedure 98 
5.3 Further Examples of System Design 105 
5.4 Sampled-Data System without a Sampler in the 
Error Channel 111 
5.5 Systems with Multiple Inputs bie 
5.6 Suggestions for Further Work 126 
VI DESCRIPTION OF THE ANALOG COMPUTER SET-UP. 128 
VII BIBLIOGRAPHY i 141 








lv 





CHAPTER I 


INTRODUCTION 


1,1 General Background 


The past few years, since the end of World War II, has witnessed 
a tremendous increase in industrial and military applications of feed- 
back control systems. Asa result of this activity, the field of feed- 
back control systems has developed and specialized in many direc- 
tions. One important class of systems that has recently been sub- 
jected to vigorous investigation is known as sampled-data systems( or 
sampling servomechanisms), In this class of systems one or more .of 
the variables consists of trains of pulses or sequences of values, These 
pulses are generally available at equally spaced intervals of time and 
no information is received between two consecutive pulses, Conven- 
tional systems, or systems in which the variables are always continu- 
ous functions of time, may, in fact, be considered as a special limit- 
ing case of sampled-data control systems, where the sampling interval, 
or period, approaches zero, However, the present body of published 
sampled-data control system theory is most useful when applied to sys- 
tems that have sampling rates low compared to the maximum response 
frequency of the system. It is to that type of system that the findings 
presented in this manuscript are directly applicable, 

One of the earliest applications of sampled-data systems was a 
tracking system including a search radar, In this case the antenna, 
rotating at a constant angular velocity, results in information about 
the location of the object which is available only as the antenna sweeps 
past the object. Time multiplexing of telemetering data is another com- 
mon source of sampled-data, The increasing trend toward systems 
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containing a digital computer as one element is furnishing further 
examples. 

In the aforementioned applications the processes require that the 
information be in sample form, Inother cases sampling is purposely 
introduced although the signal information is available as a continuous 
time variable. Generally this is done for the express purpose of 
increasing sensitivity or for more economical use of equipment. The 
electromechanical galvanometer and chopper bar is a classical exam- 
ple of such a device (1). 

Among the early contributors to the study of these systems were 
Oldenbourg and Sartorius (2). They applied difference equations to the 
analysis of sampled-data control systems. The solutions of these equa- 
tions yielded the absolute stability of the system as well as the response 
to various inputs, This approach offered the same difficulties as exper- 
ienced in the classical methods of using linear differential equations 
for describing continuous signal feedback control systems. The method 
is cumbersome and does not permit relating the performance of the 
systems to its constants, Asa result of these difficulties, attention 
has been given to developing transform calculus methods similar to 
those available for continuous systems, 

Through the independent efforts of Hurewicz (3), Stibitz and Shan- 
non (4), transform methods were evolved which gave input and output 
relations of the system as a function of z, where z is defined as tan 
(s is the Laplace complex variable.) This actually is a special case 
of the Laplace transform in which the transform is a rational function 
of z . Because of this, the inverse transform yields the time solu- 


tion only at the sampling instants, 














Ragazzine and Zadeh (1), pursuing the developments of Hurewicz, 
formalized a transform calculus now known as the z-transform, Trans- 
form pairs similar to Laplace transform pairs were made available. 
Transfer functions in terms of the variable z were shown to be deriva- 
ble in straight forward manner. Even more important, they demon- 
strated how sampled-data could be analyzed by z-transform methods, 
This article laid a firm foundation for future frequency analysis inves- 
tigation. Significant among the many recent contributors are Jury (5) 
(6) (7) (8) and Lago (9) (10). They showed that methods and concepts 
associated with the analysis of continuous systems could with some 
modification be directly applicable to sampled-data control systems 
using the z-transform, Correlation was established between the fre- 
quency response and the transient response of some second order sys- 
tems. This information was shown to be useable for approximate design 
of higher order systems where dominance holds, 

Linvill (11) showed that the transfer function in a sampled-data sys- 
tem could be represented in the frequency domain by an infinite series, 
Fortunately, only two or three terms are generally needed for a good 
approximation, The labor involved is about the same as with z- trans- 
form techniques but it has the advantage of being easily applied by 
those familiar with the standard methods of analyzing continuous vari- 
able systems. 

As the previous discussion indicated, most investigators of sam- 
pled-data control systems have concerned themselves with analysis. 
They have shown that either the z-transform or Laplace transform 
methods can be applied to the open-loop system and readily yield infor- 
mation as to the response of the closed-loop system, Although the 
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labor involved is considerable greater than for continuous systems of 
the same complexity, both methods are straight forward, That the 
methods are applied to the open-loop system is significant, since in 
general this is known in advance by the designer, Furthermore, com- 
pensating networks, employed to improve the closed-loop responses, 
are added as part of the open-loop system. To date, design procedures 
are largely based on analysis and trial and error methods, Predic- 
tion of high order system time domain performance is based on study 
of second order systems, Inthese respects the design of sampled- 
data systems is similar to the commonly employed methods for design- 
ing continuous closed-loop systems, 

In continuous system design considerable use is made of the abi- 
lity to subdivide the system into transfer function blocks, the require- 
ment:being that each block not affect the operation of the adjacent blocks, 
This requirement is satisfied by many physical devices and consequently 
the open-loop functions of continuous servo systems can readily be rep- 
resented by tandem transfer functions, An analogous situation exists 
with respect to sampled-data control systems, However, in this case 
the transfer function must be applied to all components between output 
and input of samplers, If a system employs only one sampler then one 
z-transform function must be used to describe the entire open-loop 
system, Whenever a new network, such as compensating network, not 
employing a sampler, is introduced into the system an entirely new 
transform must be recalculated for the complete system, Even the 
slightest change ina parameter value of an existing component will 
require the determination of a new z-transform, (Fortunately, the 
introduction of a sampler into an otherwise linear system does not 


affect linearity; consequently, gain enters into the z-transform asa 
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multiplying factor, ) Because of these complications, the locus shap- 

ing on the Nyquist diagram or Bode diagram requires prodigious amounts 
of labor when compared with the amount of work required for a simi- 

lar analysis of continuous systems, If this complication could be removed 
design of sampled-data control systems would be no more difficult than 
for continuous ones, 

This investigator spent considerable effort studying this problem 
in hopes of finding some design methods which would circumvent or 
ameliorate this complication, It readily became apparent that the solu- 
tion was not to be found by changes or additions to present methods 
which modify the open-loop transfer function, Further progress could 
only be had by adopting a completely new approach, 

Since working with the open-loop transfer function necessitates 
considerable labor and complications, it was thought that perhaps the 
closed-loop transfer function might offer opportunities for design with- 
out the proportional increase in labor experienced when this approach 


is used in continuous systems (12). Since system performance is gen- 
, es ieneneel 
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erally specified in terms of closed-loop performance, this approach 
would actually have the advantage of being a design method based on 
systhesis as contrasted with the present methods based on analysis. 
Investigation along this line produced a direct synthesis method which 
permits the designer to shape the time response of the closed-loop 
system. Also, the procedure indicates and helps to resolve possible 
conflict between the performance specifications and the plant (fixed 
components inthe system), When applied to higher order systems, 
this method requires less time than those methods whose starting point 
is the plant transfer function, The detailed results of this investigation 
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constitutes the contents of this dissertation. 


1,2 Discussion of the Dissertation 

After a few introductory comments the subject material of this dis- 
sertation will be discussed by chapters, The figures showing system 
time responses were determined by an analog computer, These com- 
puter curves were found to check out the analytical calculations for 
values at the sampling instants exceedingly well, the error, in general, 
being less then reading errors of the 4 centimeter paper. So far as 
possible, the nomenclature used in this dissertation is that recom- 
mended by the AIJEE Standard Subcommittee on Terminology and Nomen- 
clature of the Feedback Control Committee (13). Because of the nature 
of this dissertation additional terminology was required, 

Chapter II presents various aspects of current sampled-data theory 
which will be required for subsequent developments in this disserta- 
tion, Because of the newness of the field, important developments are 
found scattered throughout the literature and it would not be expected 
that the average reader have available all the references from which —— 
the material of this chapter is drawn, Hence, Chapter II presents this 
material in more detail than would be normally expected, For the rea- 
der familiar with the field, this chapter could be scanned or omitted 
completely. No new contribution is made in this chapter but the author's 
explanations may add further insight into the theory, 

Chapter III is the major contribution of this dissertation, A philo- 
sophy of sampled-data system design is developed which is nearly pure 
systhesis, With this method the designer is able to take time domain 
specifications and proceed ina logical manner to the final system, 
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The method permits including specific components as part of the spe- 
cifications, 

The procedure is quasi graphical leaning strongly on root locus 
techniques (14), This will find favor with practicing engineers who 
generally prefer graphical over analytical methods since they are less 
subject to mistakes and usually provide more circumspection, 

The procedure may be described by the following five steps: 

(1) The poles and zeros of the closed-loop function are deter- 

mined from the time domain specifications, 

(2) The poles of the open-loop function are found by graphical 

methods, 

(3) The open-loop poles and zeros are shifted to include the 

important poles and zeros of the specified components. 

(4) The new closed-loop poles and zeros are determined gra- 

phically and checked against the specifications, 

(5) The appropriate compensation function is determined, 

Steps (3) and (4) may be eliminated if the designer feels that 
the system will operate satisfactorily with a compensator cancelling 
out all the poles and zeros of the specified components, 

Chapter IV develops the necessary techniques to implement the 
design philosophy described in Chapter III Step (1) requires a set of 
correlation theorems which permits interpretation of the closed-loop 
system time responses, such as peak output, steady state error, set- 
tling time, etc, in terms of the pole-zero configuration of the closed- 
loop system frequency function. Some correlation theorems are avail- 
able in the literature (6) (15). However, the set of error. coefficient 
theorems and the settling instant theorem presented in this chapter 
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were derived by the author, 

The terms closed-loop root locus and open-loop root locus have 
been introduced to distinguish between the former which is the well 
known ''Evans Technique" (14) and the latter which is a method of 
finding the open-loop poles from the closed-loop function, In this dis- 
sertation the open-loop root locus was a necessary development to 
implement step (2). This technique, originated by the author, repre- 
sents an important addition to the present root locus methods for sam- 
pled-data control systems, Step (3) presents no difficulties while step 
(4) utilizes the standard root locus techniques, Current literature has 
numerous instances where the closed-loop root locus or Evans locus 
has been applied to sampled-data control systems (6) (12) (15), Only 
portions of these references which are directly applicable to this design 
method are included. The details of step (5) are left to Chapter V, 

Basic pole-zero configurations are derived which give the engin- 
eer a starting point for his design, It is shown that most specifica- 
tions can be satisfied by a basic pole-zero configuration consisting of 
one pair of complex dominant poles plus several additional sets of real 
poles and zeros which must be carefully located according to the rules 
developed in this dissertation, 

Chapter V applies the technique developed in Chapters III and 
IV to a wide variety of sampled-data control systems which included 
hold networks and delays. Although not an exhaustive list, the success 
of the design methods on these problems indicates its general applica- 
bility. The steps involved in deriving the Laplace transform from the 
z-transform are illustrated, To the author's knowledge this is the 
first time that this process has been presented in such detail. The 
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experimental findings verifying the methods of this dissertation are 
presented, 

The general theory is extended to types of sampled-data control 
systems which do not enjoy a closed-loop transfer function inde pendent 
of the input function, This is a further complication found in sampled- 
data systems not having a direct analog in continuous linear systems, 
These systems are physically characterized by the omission of a sam- 
pler in the error channel, Conventional analysis methods readily test 
the absolute stability, However, to the author's knowledge this disser- 
tation is the first to include design procedure for these very special 
types of control systems, 

Methods for handling systems with multiple inputs », suchasa 
disturbance input, are also discussed, 

Chapter VI presents the analog computer set-up that was used 
to verify the methods in this thesis, The use of an analog computer, 
with the transfer function of the various components making up the sys- 
tem directly analoged, permits the system to be analyzed without the 
effects of nonlinearity overshadowing the features of the theory to be 
tested in the experiment. A specially designed mechanical sampler 


(16) and several new experimental techniques are discussed, 





CHAPTER I 


INTRODUCTION TO SAMPLED-DATA THEORY 


2.1 Operation of the Sampler 

A linear sampled-data feedback system can be represented by the 
introduction of a switch in the otherwise continuous linear system, as 
shown in Figure 2.1. This switch, referred to as the sampler, is the 
component that distinguishes a sampled-data system from a continu- 
ous system. Figure 2,1 is misleading regarding the operation of the 
sampler, for in physical systems the sampler operates in accordance 
with the circuit shown in Figure 2,2. Here the switch, SW, is alter- 
nately closed on position a for T. seconds and closed on position b 
for Ty seconds, The travel time between positions a and b is con- 


sidered negligible and the sampling period T becomes 


T = Ad (2. 1) 


b 
In this dissertation the sampling period is assumed constant, If more 
than one sampler is present inthe system, they are considered as 
operating with the same period and synchronized, 

With the sampler illustrated in Figure 2,2 , the output signal is 
other than zero only for the T. seconds each period the switch is in 
position a , Thus the operation of the sampler converts a continuous 
signal e(t) at its input into an output signal of pulses with width qT. 
and a height proportional to the amplitude of e(t) at the sampling instant, 
Figure 2,3 illustrates the input and output signals of a physical sam- 
pler, 

In the systems to be considered in this dissertation, the time T. 
is much smaller than both Ty and the time constants of the systems, 
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r(t) + 


R(s) 





E(s) E*(s) 


H(s) 


Figure 2.1, A sampled-data system 
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Figure 2,2, Physical circuit for sampler 
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O time 


(a) Continuous signal applied to the input of the 
sampler 


time 
O 
T 
(b) Pulse train at output of the sampler 


Figure 2,3. Input and output signals of a physical sampler 














Physically this latter requirement means that the system responds a 
negligibly small amount while the pulse is applied, If this is the case, 
the system response would be the same to an impulse with strength 
equal to the area of the pulse. Since mathematically the impulse is 
much simpler to handle than the pulse, both in the literature and in 
this dissertation, the impulse train is substituted for the actual value 
of the continuous time variable, It should be noted, that the area of 
the impulse is considered as equal to the value of the continuous vari- 
able at the instant of sampling rather than corresponding to the area 
of the pulse. With the area of the pulse equal to the value of the con- 
tinuous variable multiplied by %. , the mathematical representation 
of the sampler inveltes a gain change, Thus, the gain K of the open- 
loop system is related to the actual gain Ko used in the physical sys- 


tem by: 
K, es (2, 2) 


The Figure 2,4 may yield more insight to the nature of the approx- 
imation involved when the pulse is replaced by an impulse under the con- 
ditions stated above, It is clearly shown in this figure that the fre- 
quency content of the pulse is nearly constant over the extent of the 
bandwidth of the open-loop function, Thus, the system output contains 
essentially the same proportion of frequency components with either 
the pulse or impulse input, The change in gain is illustrated with T. 
less than one second, 

It is helpful for analysis to consider the sampler as an impulse 
generator which is modulated by the input signal, The process of con- 
sidering sampling as modulation of an impulse train is shown in Fig- 
ure 2,5 , wherein the values of the impulses are represented by the 
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(a) The pulse and impulse time functions 
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(b) The frequency spectrums of the impulse (1) , pulse 
(2) and bandwidth of the open-loop system (3) 


Figure 2,4, The impulse approximation illustrated in terms 
of the Frequency Domain, 
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(a) Input to sampler 


i(t) 
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(b) Impulse train 
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(c) Output of sampler 


Figure 2.5 . The sampler as an impulse modulator 





the amplitude of the lines occurring at the sampling instants, With > 
input and output denoted as e(t) and e*(t) respectively, the output 
may be written as 

e*(t) = e(t) i(t) (2,3) 


where i(t) represents a train of unit impulses 


OS 
i(t) = y u(t - nT) (2. 4) 
=—- O© 


Equation 2,3 may be written equivalently as 


oO 


e*(t) = oll e(nT) Uo (t - nT) (2, 5) 
n=o 
Here the negative values of n are omitted by virtue of the assump- 
tion that e(t) is zero for all negative values of time, The Laplace 


transformation of both sides of Equation 2,5 yields 

CO 

- nTs 
E*(s) = ) e(nT)e (2. 6) 
n=o 
An alternate and very useful expression for E*(s) can be 

obtained by the use of contour integration, Equation 2,3 represents © 
the sampled time function as the product of two time functions from 
which the Laplace transform can be obtained by a process of complex 


convolution, The complex convolution integral is given by (17) 


c+je 
L[f, (t) f, (t) ] = zr | FY (s - w) F, (w) dw (2, 7) 
c- joo 
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The contour of the integration of Equation 2,7 is a pathfrom c - joo 
to c + joc along an analytic strip between the singularities of 
Fy (s - w) and those of F, (w) . 


Letting the impulse train be represented by fy (t) of Equation 


Ze, Fr (s - w) becomes 
oo 
F(s - w) = I(s - w) = \ e- nets) (2, 8) 
n=o 


5 


This infinite summation can be expressed in closed form as 


ie - ee (2.9) 


1 - e7 T(s-w) 


Then, equating e(t) to f, (t) yields 


F,(w) = EF (w) (2, 10) 


With the indicated substitutions in Equation 2,7 the Laplace transform 


of e*(t) becomes 


ct#+jco 
L [e* (t) ] = ar | E (w) Te om (2, 11) 
= S-W 
€ - joo 


Evaluation of Equation 2,11 can be effected by closing the poles 
I (s-w) in the right half of the w-plane as shown inFigure 2,6, When 
the degree of the denominator of E(s) in s is higher than that of the 
numerator by at least two, the integral is zero along the infinite semi- 
circle, The alternate path in the left hand plane will be considered in 
a later section of this chapter. 


Applying Cauchy's integral formula gives 
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C tjJQ0 w-plane 


Poles Of Ew) 
yy 


C-j® 


Figure 2,6, Paths of integration used to evaluate the 
complex-convolution integral, 





18 





3 
rm (<a dw = » Res E(w) (2212) 


l-e = T(5-w) poles of aver DAEEas) 


I(s-w) 


The poles of the I(is - w) are the zeros of 





le "ee UNE (2. 13) 
Or 
Le/ Bog OMe 2urn (2. 14) 
n T 


These poles, though infinite in number, are all simple and hence the 


residues are given by Expression 2,15 
n= oo 


” eres) re (2, 15) 
Ee ices il oa Soe 


Since 


21m 


djaw Fe 2 ~Tel 2 Tn (2. 16) 
O9= 8+ J—y— 


The value of the integral represented by Equation 2,12 is 





oO 
E*(s) = =) E(s+j—2*) (2,17) 
n=- oo 
Defining the sampling frequency as Le whe re 
Gy = 2n/T (2, 18) 


Equation 2,17 may be written as 
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n= &© 


ca << —— E(s + jnUJ,) (2. 19)# 
T 


N=-o0 


It is seen that two equivalent expressions are available for analy- 
tical description of the sampler, as demonstrated by Equation 2,16 
and Equation 2,19, The first equation, which expresses the Laplace 
transform of the sampler output in terms of the values of the input 
function at the sampling instants, is particularly helpful in gaining 
insight into the physical operation of the sampler and its effect on sys- 
tem performance, The second equation places in evidence the fre- 
quency effect of the sampler on the input signal, 

Replacing s by s+ j mis. , where m is any integer, in Equa- 
tion 2,19 or its alternate expression Equation 2,6, results in an 
expression which is identical with E* (s) (1) (18) . Thatis, E* (s) 
is a periodic function of s with period J Ud, . Expressed symbo- 
lically ; 


E* (s + jml/,) = E* (s) (2, 20) 


This indicates that the s-plane can be divided into strips of equal width 
with centers at n(J, as shown in Figure 2,7, 

The spectrum alteration performed by the sampler is shown in 
Figure 2,8, From this it is seen that the sampler contains an exact 


replica of the input as well as an infinite number of additional images 


*Equation 2,19 is in error if e(t) # 0 at t= 0+ . In this case the 


Equation 2,19 should be modified to read . 
n= co 
l 
E* (s) = - \ (s + jn L/ ,) + 1/2 e(0+) 
n= - oo 
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Figure 2.7. Period strips of E* (s) 
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(a) 


(b) 


Figure 2, 8, 
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Amplitude spectrum of E(s) 
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Amplitude spectrum of E* (s) 


Input and output amplitude spectrum of 
a sampler, 
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centered about a multiple of the sampling frequency. These additional 
components are referred to as complementary signals (11). Spaces 
between images would indicate that it is possible to regain the exact 
signal by filtering, Unfortunately, a satisfactory filter for the process 
would introduce too much time delay. Any Ee with frequencies 
greater than the sampling frequency would overlap into the adjacent 
strip with the result that the information would be lost, This conclu- 
sion is equivalent to Shannon's sampling theorem which requires that 


the sampling frequency be at least twice the highest input frequency. 


2.2 Elements of z-Transform Theory 

If the path of integration for Equation 2,11 had been closed to 
the left, the finite number of poles of E(w) would have been enclosed 
and there would be a finite number of residues to evaluate. This alter- 


nate path is also indicated in Figure 2,6, In this case 


L [e* (t)] = 3 Res E(w ) ee (2,21) 
~T(s-w_) : 
poles of l-e = 
E(w) 


The use of Equation 2,21 may be illustrated by determining the 
z-transform of the sampled sequence of a unit ramp, 


e(t) = ¢ (Z, 22} 


4 s @ 


The addition is a result of the difference in definition of the sam- 
pling process in Equation 2,6 which implies a full sampling at t=0 
and Equation 2,11 which assumes the impules to be an even function 
sampling only one-half of the value e(o). In feedback control sys- 
tems Equation 2,19 is normally adequate, 
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The Laplace transform is then 


Bis) ve a (2, 23) 


S 
E(s) is seen to have a second order pole at origin, Performing the 


operations indicated by Equation 2,21 yields 


sT 
E*(s) = ie ; 
sT am 


(2, 24) 





(e 
An alternate method of deriving Equation 2,24 employs Equation 2.6 


asa starting point. Hence 
e*(t) = Tu(t - T) + 2Tu(t - 2T) + 3Tud(t - 3T) +... (2,25) 


The Laplace transform of this series of samples is 

















E(s) = em + Maye + Te ' ar ee (2, 26) 
Multiplication of both sides of the above equation by yields 
E*(s) . l + I + Shire ele (2,27) 
Ts 2Ts 3Ts 4Ts 
Te e e e 
Multiplication of both sides of the above equation by dem and 
integration gives 
x 
Eee e@apee faepe | os ge ee, 8S (25:28) 
sT Ts 2Ts 3Ts 
Te e e e 


Where C_ is the constant of integration, The series on the right of 


the above equation may be rewritten in closed form, 
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pBisee 4.8) gy ee eee (2, 29) 
T eat esl 24 
Differentiation of Equation 2,29 with respect to Pe. yields 
8 T 
Bee) tl ee (2. 30 
eet ‘ 1)7 
which is identical to Equation 2,24, The substitution z = Md into 


the expressions gives rise to the z-transform, The above expression 


in z-transform notation becomes 


a (2. 31) 


(2 - 1)" 
3 
By these methods it is possible to derive a table of z-transforms 
(5). Table 2,1 is a short list of z-transforms, 
The inversion of a z-transform function to obtain the time res- 
ponse of the function at sampling instants can be effected by anyone 


of three ways: 


(1) Evaluation of the inversion integral (19) 


c(nt) = a Cine ae (2, 32) 
[~ 


where [° is a path of integration enclosing all the singular 
points of C(z) 

(2) Partial fraction expansion of the function so that each term 
can be identified with the corresponding time function of 
Table 2,1, 


(3) Series expansion in descending powers of z in which the 
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Table 2.1. Short List of z-Transforms 


Tinve fun cTien Lap/ace Z-lranstorm 
Transform 





* aomay be complex 
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coefficient of 1/z” corresponds to the value of the time 


function at the nth sampling instant, 


2.3 z-Transfer Function 
Of primary importance in servo analysis are the input-output 
relations of networks, This relation for sampled data systems will 
now be established using the simple systems illustrated in Figure 2,9, 
which shows a linear network W _ pulsed by a sampler, In addition 
to the desired continuous output time signal a fictitious sampler is pro- 
vided to generate a sampled output signal, Although not actually part 
of the signal flow in the system, consideration of the sampled output, 
rather than the continuous output, simplified the analysis as will be 
shown, 
It is evident from Figure 2,9 that the Laplace transform is 
given by 


C(z) 2 G(s) E* (s) (2, 33 


where G(s) is the Laplace transform of network W and E* (s) 
is the sampled input. Applying Equation 2.19 yields the transform 


of the pulse output of the output sampler. 
= l b ote ° 
C#(s) = a) Gls +in Wj) E*(s+jn WJ.) (2.34) 


Recalling the periodicity of pulsed-transfer function as expressed by 


Equation 2,20 , the above expression may be rewritten 
l : 
ae = % 
C* (s) a E* (s) G(s + jn cy) (2, 35) 
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Figure 2,10, A system witha sampler between 
networks, 
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Since E*(s) is a factor in each term of the above series, it may be 
placed before the summation, hence 
| oO 


Cx (s) = E* (s) lar) le +n Wy, (2, 36) 


n= - oo 


Denoting the bracketed term by G* (s) 


oO 
G*(s) = ae y G(s + jn U/) (2, 37) 
CO 


n=- 
Equation 2,36 may be rewritten 
C#(s) = E* (s) G* (s) (2, 38) 
In the z-transform notation 
C(z) = E(z) G(z) (2, 39) 


It is seen that the z-transform of the output of a network is simply the 
z-transform of the input multiplied by the z-transfer function, Fur- 
thermore, the z-transfer function is derived in the same manner as 
the transform of a time function. 

An alternate expression for G*(s) interms of g(t) , the 
impulse response of W , can be derived by noting the similarity in 
form between Equation 2.19 and Equation 2,37 given below for com- 
parison 


oo 
E*(s) = yy E(s + jn UJ.) 


nse oS 


oO 
G*(s) = ae y. Gie + jn LJ.) 


nz=- OO 
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Taking the Laplace transform of Equation 2.3 yields 


E* (s) = L[e(t) x i(t) ] (2, 40) 


Since G(s) is the Laplace transform of g(t) it follows at once that 


G* (s) = Lf g(t) * i(t) ] (2, 41) 


The over-all transfer function of two networks separated by sam- 
plers is readily seen to be equal to the product of the z-transform for 
each network, This is illustrated in Figure 2,10 , If no sampling 
device separates the networks, the input to the second network is the 
continuous function of time present at the output of the first network, 


Then the corresponding pulsed (starred) transfer function will be 
| nz ©° 


G* (s) = ny G{s + jn / ) G,(s +jn UJ) (2, 42) 


n=z=-0Co 
Hence the over-all transfer function must be calculated by considering 


the system as a whole, The accepted notation (1) for this case is 


G* (s) = G, G*, (s) (2, 43) 


or in terms of z-transform notation 
G(z) = G, Gy {z) (2, 44) 


The transfer function of a system with feedback, such as shown 
in Figure 2,11 may be derived in the following manner, The expres- 


sion for the error is 
E(s) = R(s) -H(s)C(s) (2, 45) 
The output is directly related to the sampled error by 
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Figure 2,11. A simple feedback sampled-data system 
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C(s) = E* (s) G(s) (2, 46) 


Equation 2.46 substituted into Equation 2,45 yields 
E(s) = R(s) -H(s)G(s)E* (s) (2,47) 


Taking the starred transfer function of both sides 


Ex (s) = R* (s) -HG* (s)E* (s) (2, 48) 
Solving for E* (s) 
Et (sje, —_ el) 2 (2. 49) 
1 + HG* (s) 
The output is then 
C* (s) = ___R* (s) G(s) (2, 50) 
1 + HG* (s) 
Taking the starred transfer function of both sides 
Geely, Se eats! (2, 51) 
1 + HG* (s) 
and in z-transform notation 
C(z) = __R(z) G(z) (2, 52) 


1 + HG(z) 


Table 2,2 gives the transforms of the output for a number of 
basic servo configurations. Table entries 2, 6 and 7 show the input 
may not always appear as a simple multiplicative factor in the expres- 
sion for C(z). Inthis case it is not always possible to derive an 
expression for the feedback system transfer function independent of 
the input. As discussed in the introduction this creates a type of sys- 


tem whose analog is not found in linear systems operating wholely on 
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Table 2.2, Output Transforms for Basic Sampled-Data Systems 
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continuous data, 


2,4 Final and Initial Value Theorems, Error Coefficients and Mapping 


The output of a system C(z) may be expanded in partial fractions 





as 
Dz Dyz 
C{z) * —— $+ ype tt ay tee (2. 53) 
Zz-l as | "ie 
Z-e Z2-€ 


The inverse transform into the time domain is readily determined by 


entry 4 in Table 2,1 


-d,T -d,T 
c(t) = D + D 


a if Scacice (2, 54) 


Assuming that the real parts of the d's are positive all terms of Equa- 
tion 2,54 goto zero as time goes to infinity with the exception of D, P 
Therefore, the final value of c(t) may be found directly from C(z) 


as follows 


lim(z - 1) C(z) = lim c(t) (2, 55) 
zZ—>!l t —>cc 


Expanding C(z) in descending powers of z gives 


nal H. 
C(z) = H + + =o ae (2, 56) 
Z 





The inverse transform into the time domain is 


c(t) = H, + H,Uj(t-T) + H,U(t ~ 2T) +... (2, 57) 


Hence, Ho , the initial value of c(t) , may be found directly from 


C(z) by 
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lim C(z) = lim c(t) (2, 58) 
Z— oo t >0 


The error expression for single loop systems with unity feed- 


back is 


ini eee oe (2, 59) 
1 + G(z) 


For a step input 





Riz) = ——— (2, 60) 
z-il 
Thus 
E(z) = —————_____ (2, 61) 
(z-1) [1 + G(z) ] 
Applying the final value theorem the steady state error eae is 
e.. = lim 22) | —__-2_____ | = lim —_—__(2, 62) 
hee (z-1) (1 + G(z) ) Ze Gz) 


Then the error coefficient for a sampled-data system may be defined 


as 
Kx = Controlled Variable (2. 63) 
P Actuating Error 
or 
l-e 
K* ee ee Se lim G(z) (2, 64) 
P Or z—>l 


In like manner it can be shown 


Kx = ee linn te eee (2, 65) 
zl 
and 
i ee limeizenic Gian (2, 66) 
c — 
i z—>l 


It was explained in Section 2,1 that the s-plane could be consid- 
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ered as divided into strips, The conformal mapping of the first strip 
in the s-plane is the unit circle in the z-plane, This and mapping of 


other important lines are shown in Figure 2.12. 


2.5 Response between the Sampling Instants 


The output of a sampled-data system is continuous even though 
its input is sampled, but the z-transform only contains information of 
the time response at the sampling instants, This section presents an 
extension of the z-transform called the modified z-transform which 
permits determining the output of a network between sampling instants. 

Consider the system in Figure 2,13 which introduces a ficti- 
tious delay between the output of the network and the sampler at the 
output, A represents the fraction part of the delay to the total sam- 
pling period and hence has a positive value between zero and one, The 
effect of this delay on the signal presented to the sampler is illustrated 
in Figure 2,14. For reasons which will become clear later it is aus- 
picious to define the change of variable 

OS = l-m (2, 67) 
E juation 2,67 implies that mT is a pure advance. A system incor- 
porating this pure advance before the sampler and a delay of one sam- 
pling period following the sampler, as shown in Figure 2,15, is seen 
to have the same output as the system in Figure 2,13. The impulse 
response of the network and pure advance is g(t + mT) where a(t) 
is the impulse response of the network alone, Employing Equation 


2.41 it follows: 


Cy (z) = L[ g(t +mT) x i(t) ] (2. 68) 
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Figure 2,13. Sampled-data system with fictitious delay, AT. 
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Figure 2,14, The output compared with the delay output of 
system in Figure 2,11. 
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Figure 2,15. Sampled-data system with output equivalent to 
that of system in Figure 2,13 , 


og 


Equation 2,68 may be evaluated using the complex convolution integral 


in the following form: 


“Ctjcs 
Co (z) = wae G(w) ut meee ee dw (2, 69) 
c-j Co 1 - o7 oe 


Evaluating the integral in the left hand plane provides a closed 
form for C* (z) . It should be noted that the substitution represented 
by Equation 2,67 assured the convergence of the integral on the infin- 
ite semi-circle, 

Because output C, (t) is identical with C* (t) but shifted in 


time by one period, it follows from the real translation theorem (17) 
-sT 
C* (s) = e C* (s) (2. 70) 


Varying the value of m from zero to unity, the entire response 
between the sampling instants can be found, The z-transform evolved 
from this technique is called the modified z-transform, With this nota- 


tion the transfer function of the system in Figure 2,13 , including the 


linear network and delay AT, is 
ms T 
G(z,m) = 2! yy Residue | —oi8) (2, 71) 
poles of CS ae 
G(s) 


Table 2,3 is a short list of modified z-transforms for various 


forms of G(s). 


2.6 Analysis and Design of Sampled- Data Control System by 


Conventional Techniques 


The absolute stability of sampled-data control systems can be 
40 





Table 2,3. Short List of Modified z-Transforms 








Laplace Transform Modified z- Transform 
/ Pee 
/ oa 2 
Le mT eZ 
ro sa 2 -/ (2 -/)% 
=O1E 
Kia / Cc 
S+q 2- CW 


* a may be complex 


41 





determined in much the same manner as is done with continuous sys- 
tems, For continuous systems a very powerful technique has been 
devised involving the construction of asymptotic diagrams in which 
log modulus and phase are plotted against log frequency. This is 
known as the Bode diagram and is primarily useful when the transfer 
function is expressed as a product of simple factors and the variable 
increases indefinitely in magnitude. That latter requirement is not 
characteristic of the z-transform, where for this diagram, z is 
varied along the unit circle. 

A Nyquist diagram may be constructed by plotting the gain and 
phase of the open-loop transfer function as z is varied along the unit 
circle, Closure of the coordinate -1+ j0 indicated instability. Con- 
struction of appropriate oe circles (20) gives an indication of the 
relative stability. 

As an example of analysis and design using the Nyquist diagram, 
consider the system shown in Figure 2,16, The sampling period is 
0.5 seconds, The first step in the analysis involves determining the 
open-loop z-transfer function, G(z). This requires a partial fraction 
expansion of G(s) , converting each term to z-transform with the 
aid of Table 2,1 and finally a collection of these terms over a common 


denominator, These steps are illustrated below, 


10 





G(s) = ——_—— (2,72) 
s(s + 3,2) 
G(s) ae. (2. 73) 
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a a G(s)=_ —/2 


ze S( 543.2) 


Digure 2,16, Sampled-data system to be compensated. 
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G(z) = pp gy (2, 75) 


Analysis on the Nyquist plot, Figure 2.17 shows the system to 
be unstable. Let the design require a peak overshoot of 50% in the 
time response while maintaining the same low frequency gain, Asa 
first try, it will be assumed that this requirement could be met with 
phase lag network compensation 

2 st+Z 
Sop MO aS ores an Og oath 
The new z-transform open-loop transfer function must be deter- 


mined as shown in Equation 2, 77 through 2,80. 


G(s) = See FST —_ 
G(s) = =e. BA Se C278) 
G, (a) ¢ Se1288e . 18s | 16552 (2 9) 
G, (z) = a (2, 80) 


The resulting Nyquist diagram, Figure 2,17 , gives the locus 
touching the Soe = 2,5 circle, This indicates a time response over- 


shoot of 53% . Trying a second network compensator 


Gg) 18544 2.81 


The G(s) and G(z) of the new open-loop system are then 


G(s) = qee aeas (2. 8%) 
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G, (z) = oe Sey (2, 83) 


Z — Ze e Z2ae- e 


The resulting Nyquist plot, Figure 2,17 , shows the locus to be 
outside the Mes equal two circle, indicating a peak overshoot of less 
than 44% . Hence, the design is completed as far as this procedure 
is concerned, The first step in calculation of the time response at 
sampling instants is to form the output z-transforms as given in entry 
4 of Table 2,2, This is done in two steps, first by forming the closed 
loop transfer function and then followed by multiplication of the trans- 


form of the input, The final result is 


2 
C(z) = ven eb SEM (2. 84) 


oe = ola + 682° ~ ,46z + .121 


The first few values of the time response at the sampling instants 
are placed in evidence by dividing the denominator into the numerator 


as 


lao 7 + 1,35 . 76 1, 03 


C(z) = ey = a Fees (2, 85) 





Z 


Although the actual peak overshoot is outside of specifications, the 
Nyquist procedure has given a system which at least has a reasonably 
good time response, 

The previous design problem has illustrated the definite trial and 
error nature of design of sampled-data control systems by conventional 
techniques, For higher order systems, the labor may be prohibitive 
for designing to fixed specifications. | 

Root locus methods starting with open-loop transfer functions 
suffer from the same difficulties in determining the effect of a com- 
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pensation network, 


Examination of the open-loop transfer functions given by Equa- 
tions (2, 75) (2, 80) (2,83) permits drawing four important conclusions 
about z-transfer functions derived from continuous components, 

(1) The s-plane poles have direct images inthe z-plane 

(2) The s-plane zeros do not have direct images in the z-plane 

(3) Additional z-plane zeros are introduced when taking the 

z-transform 

(4) Ina z-transform the numerator is of one degree less than 

the denominator if the function has a zero value at zero 


time and a nonzero value at the first sampling instant, 
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Figure 3.1. A compensated sampled-data control system 
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Figure 3,2, A sampled-data control system employing a 
sampler with the compensator 
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Figure 3.1, A compensated sampled-data control system 
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Figure 3.2, A sampled-data control system employing a 
sampler with the compensator 
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culties the classical approach has been to work directly with the open- 
loop function, A compensator is chosen to reshape the Bode diagram 
according to certain tenuous relationships between the frequency and 
time domain, For continuous control systems the uncertainty in the 
final performance of the closed-loop system is compensated for by the 
ease with which the method can be applied, 

The design example in Section 2.6 showed that this ease of pro- 
cedure for design based on the open-loop function does not apply to 
sampled systems unless a sampler is included between the compensa- 
tor and the plant, Such a case is shown in Figure 3,2, It is these 
difficulties which lead to the philosophy of a design procedure with the 
closed-loop functions as a starting point. The utilization of this philo- 
sophy requires two new developments: 

(1) Method of determining the closed-loop transfer function 

from the specifications 

(2) Method of determining the corresponding open-loop trans- 

fer function, 
This dissertation evolves these two developments making possible a 
logical design procedure enabling the designer to take the specifica- 


tions and proceed ina direct path to the final system. 


iA Synthesis through the Closed- Loop Poles-Zeros Configuration 
The details of the synthesis procedure are shown by the flow 


chart of Figure 3,3, The various blocks represent steps in the pro- 
cedure, The letters beside each arrow identify the operation 
required to advance to the next block, 

The chart is headed with the time domain specifications, By 
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use of ''correlation theorems" it is possible to interpret these speci- 
fications in terms of closed-loop pole-zero configuration on the z-plane. 
From the pole-zero configuration the ''open loop root locus" yields the 
corresponding open-loop pole-zero configuration. At this juncture the 
designer can choose either of two concluding procedures, 
(1) Proceed directly to determination of the compensator 
(2) Attempt to conciliate the differences between the poles 
and zeros of the open-loop function and those of the 
plant so that at least the major poles and zeros of the 
plant are included in the final open-loop function. 
If the first method is chosen the procedure is straight forward. 
The desired open-loop z-transform is converted into its equivalent 
Laplace transfer function, The transfer function of the compensator 


may now be found by simple division 


G, (s) 
G (s) = me (3511) 


: G,, (s) 


where Ga (s) is the open-loop or "desired" transfer function cor- 
responding to the selected closed-loop function, The actual compen- 
sating network may be determined by standard network synthesis 
procedure (12). 

A study of Equation 3.1 reveals that the compensator will, in 
general, cancel out each pole and zero of the plant as well as add the 
poles and zeros of the desired function, In practice it is impossible 
to exactly cancel the poles and zeros of physical components such as 
electric and hydraulic motors since their precise characteristics are 


never known over the full range of operation, Fortunately, complete 
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Figure 3.3. Synthesis Flow Chart 
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cancellation is not necessary. System performance is not significantly 
changed if the pole and its cancelling zero only lie near each other. 
This situation gives rige to a pole and zero placed close together in the 
closed-loop function and the resulting residue of the transfer function 
in the pole is negligible. 

A more serious objection rests on the argument that it is not pos- 
sible to cancel out equipment time constants and expect the equipment 
to function in a linear mode over a large range of input values, It is 
unreasonable to expect that cancelling out a motor time constant and 
replacing it with a much smaller one will result in any increase in the 
actual maximum obtainable acceleration or that even the motor and 
gear combination will be able to follow a ramp which exceeds its max- 
imum speed, Hence for this reason it is strongly urged that the 
designer proceed according to method (2). 

In method (2) the designer compares the poles and zeros of the 
desired open-loop with that of the plant. By shifting the location of 
the closed-loop poles-zeros it is usually possible to include the major 
poles and zeros of the plant into the desired open-loop transfer func- 
tion, The open and closed-loop root locus furnish a graphical method 
for bringing about this correspondence, The impossibility of incor- 
porating the major components of the plant into the desired open-loop 
function does not imply failure of this design method, Actually, this 
is a signal to the designer that the fixed plant is not capable of pro- 
viding the necessary performance demanded by the specifications, 

This simply means that the plant must be improved or, if this is not 
practical, the specifications must be changed to concede the realities 


of the situation. Thus this method increases the designer's circum- 
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s pection, 

If the compensator can be cascaded to the plant through a sam- 
pler as shown in Figure 3,2 , the preceding steps are still taken and 
the compensator can be found directly from the z-transform of the 
system, 

Gy (z) 


G. (z) = (3.2) 
G_(z) 
p 





3.3 Time Response and Basic Pole-Zero Configuration 


The critical step in this procedure involves choosing a pole- 
zero configuration to satisfy the time requirements, The typical 
time requirements for a transient response to a step is shown in Fig- 
ure 3.4. The steady state error to a ramp and parabolic input are 
shown in Figure 3.5. Numerous pole-zero arrangements could be 
found to satisfy a particular set of time specifications, However, the 
correlations theorems to be formulated in the next chapter are based 
ona dominant pair of complex poles which together with a real zero 
will fairly well control the transient response, A pole-zero spaced 
close together and near the point 1,0 inthe z-plane can provide 
control over the steady state-error. Withthis, the basic pole-zero 
configuration in the z-plane is that shown in Figure 3,6. 

It is not to be inferred from Figure 3,6 that all specifications 
can be discharged by third order systems. However, in many cases 
the synthesis can be handled initially in this way, Asa last stepin 
the design additional poles or zeros may have to be added if the com- 


pensator is to be physically realized, (This is discussed in more 
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F igure 3.4, Typical transient response to a unit input 
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Figure 3,5, Steady state errors to ramp or parabolic inputs 
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Figure 3,6, Basic closed-loop pole-zero configuration 
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Figure 3,7, Typical closed-loop pole-zero configuration 
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detailin Chapter V.) If the poles and zeros have large real parts in 
the s-plane they will give rise to z-plane poles and zeros near the ori- 
gin and have little effect on the previously determined results, Figure 
3.7 shows atypical pole-zero configuration, 

It should be kept in mind that this control of time response 
extends only to the sampling instants. This generally is satisfactory, 
However, a complete description of the output between sampling instants 
may be secured by employing the modified z-transform technique to the 
final system, 

Plants with poles and zeros in the right hand half of the s-plane 
require special consideration, The poles will have images outside the 
unit circle in the z-plane and so will the zeros in most cases, Rather 
than attempt to cancel these poles and zeros they should be included 


as part of the final open-loop transfer function, 
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CHAPTER IV 


SYNTHESIS TECHNIQUES 


4.1 Introduction 

This chapter develops techniques to implement the design phil- 
osophy described in Chapter III], First, simple relations between the 
pole-zero configuration and the time response are developed, The 
desired response, as shown in Figure 3,4 and Figure 3,5 , may be 
specified by the following five quantities: 

(1) The peak sampling instant, aS , the sampling instant with 
the maximum overshoot for a step input, 

(2) The peak sampling overshoot, M , the amount of the 
first overshoot at the sampling instant — for a step input. 

(3) The settling sampling instant, no» the sampling instant 
that the value of the output reaches within 5% of its final 
value for a step input, 

(4) The number of oscillation, N , for the interval up to the 
Serine instant for a step input, 

(5) The steady-state error constant an : K* » and K* 
for a unit step, unit ramp and unit parabolic input respec- 
tively. 

The first four quantities describe the transient response while the last 
gives the steady-state error. 

Basic pole-zero configurations such as shown in Figure 3.7 are 

next derived giving the engineer a starting point for his design, Fin- 
ally rules for applying the open-loop and the closed-loop root locus 


are presented, 
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4,2 Correlation Theorems 
For the sampled-data feedback control system shown in Figure 


4,1 the closed-loop system function W(z) is 


Wiz) = pel = « Sh (4, 1) 

where 
Q(z) = TT (z - q;) (4, 2) 

k=] 

mm 
Pia) = TT (e-py) (4, 3) 


The values of oe and Pp; may be either real or complex, It 
is assumed that the expressions do not contain multiple roots and that 
the order of P(z) is greater thanthe order Q(z). These two 
assumptions are true for most physical systems, When the input is 


a unit step, the output at the sampling instants becomes 


c(nTy =, 271 | wee ey | (4, 4) 


= + a 


where po denotes the inverse transform, Thus employing 


Equation 2,32 


l K Q(z n 
c(nT) = aes | te zdaz (4, 5) 
where i is the contour of the integration in the z-plane that encloses 
all the singularities of the integrand in Equation 4,5, Evaluation of 


the contour integral gives 
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Figure 4.1, A sampled-data control system 


Ns 
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Figure 4,2, Important characteristics of the output 
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where 


Zz = u_+jv 
O O fe) 
d P(z) 
! — 
I= (z) a a (4, 11) 


Z=Zo 


It is important to realize that the assumption in Equation 4,10 
does not neglect all other roots but neglects all other modes because 
these modes die out quickly, The accuracy of this assumption will be 
discussed later. 

In general, the peak of the output will not occur at a sampling 
instant, If n is permitted to be noninteger, thena value of n could 
be found, so that at one sampling period of time later the output will 
be repeated, This condition is illustrated by Figure 2.2 and can be 


expressed as 
c(n +1) - ec(n) = O (4,12) 


The actual value of on is the upper integral value of the actual 
n obtained. In the derivation which follows it will be assumed that 
the solution of Equation 4,12 results in an integer value for n , 


Substituting Equation 4,12 into 4,8 and solving for 2 yields* 


Fa 


a = 7 5 - ang Q(z) + ang P(e) | (4. 13) 


* The intermediate steps leading to Equation 4,13 are quite involved, 
The more important steps are as follows: the constant term, the 
coefficient and the nth power of Z| cancel yielding 
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Then the peak sampling instant can be expressed as 


0 on = —— - (sum of angles from the zero to 
P 2 the predominate pole) - (sum of 
angles from the other poles to 
the predominate poles) (4, 14) 


cos c 0. + ang Q(z.) - ang (z - 1) - ang ae.) = 
Iz. | cos (nt+1)Q. + ang Q(z 0) - ang (zo - 1) - ang PH.) 


Expanding the right hand side of the above equation by the sum of 
angles trigonometric identity, 


cos = 0. + ang Q(z.) - ang (z - 1) - ang Pi(z.) | = 
lz | cos |» 0, + ang Q(z.) - ang (z - 1) - ang Plz) cos ‘| 


- sin n 0. + ang Q(z.) - ang (Zz. - 1) - ang PHe,)| sin | 


The following trigonometric relationships exist for the values above 





cos g = : 
lz | 
u - i 
cot fang te, - 0) | = : 
sin @ a vo “o 
: Iz _| 


Substituting these into the preceding equation yields 
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Thus, it may be concluded that additional poles increase ny and that 


a zero decreases Bs 5 
The first overshoot is obtained by substituting Equation 4, 13 


into 4,10 . Hence 


K Q(z_) | 
Moe 2 ee |z ["P sin (ane 29°) (4.15) 
Zeer 2 
[(a=t) +vee] ize) 
recognizing 
Vv 
Oo 
ake | Bae (ei |- Pe oe re (4, 16) 
NA 
Z Z 
i) tv. ] 
Equation 4,15 simplifies to 
2 we K Q(z.) F 
M = ————_—____ Ss ie |e (4, 17) 
2 2 - 
t 
[(u,-1) eye ] 1S: (z 0) 


The accuracy of the assumption that all modes other than that 
due to the constant term and the predominate poles may be explained 
readily, From the complete solution, Equation 4,8, considera cer- 


tain mode mie . This mode will decrease toa small percent of 


Tan n 0. + ang Q(z.) - ang P'(z,) - ang (z - 1) = 


Tan | + - wn, </it | 


From which Equation 4,13 follows immediately. 
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its initial value when Sy = 3 , and the coefficient D is small if it 
is not too close to the other roots, Inthe event that a certain mode 
must be accounted for the contribution of this mode at a can be 
evaluated Boys the general Equation 4, 8, 

Unless the actual value of n from Equation 4,13 is equal to an 
integer an error is introduced giving an optimistic value of M. This 
is because the value of the sine term in Equation 4,15, for the non- 
integer n , is less than the sine term would be at oe . his error 
is small unless the sampling period is a large proportion of the period 
of oscillation and n is near midway between Peointeten values, 

For a zero displacement error system the constant term must 


equal unity when the system responds to a unit step of position. Hence 


m 
ai 
i=] (1 = P;) 
1 (ey, Nees ake (4, 18) 
Ww 


Ta- 
et %) 


Substituting Equation 4,18 into 4,17 the peak sampled overshoot can 
be expressed as 


M = (Product of distances from all poles to the point 
of 1,0), excluding distances from two predomin- 
ate poles to the point 1.0)/(product of distances 
from all poles to predominate pole u_ + j 9 
excluding the distance between predominate 
poles) x (product of distances from zeros to the 
predominate pole u_ + jv_ )/(product of dis- 
tance from all zeros to the point 1, 0) Zo "Dp 


(4,19) 
It is evident from Equation 4,19 that the smaller ae the 


larger M, Thus the choice of Dy and M requires a compromise, 


The effect of real poles on M_ is evident from Figure 4.3, For 
bie 
position one near the origin = is nearly unity and the effect 
l 
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Figure 4.3, The effect of a real pole on M and a 





Figure 4.4. The effect ofa real zeroon M_ and 
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oe are 
e ° e ¢ Z 
of the pole is slight. However, at position two, nn | 
will increase. Finally, a pole in position three near c , 


b,c 


b3 a 
clusions may be deduced when the poles are complex. 


<) &) M is reduced while a is increased, Similar con- 


The effect of a real zero at various positions is evident from 
Figure 4,4, As the zero approaches the point 1,0 the ratio becomes 
larger and —— increases while reducing we . A zero placed 
very near the point 1.0 is not desirable since M _ increases tre- 
mendously. For most cases, it has been found that one zero placed 
between the origin and a position somewhat less than the real part of 
the predominate pole will result in a satisfactory transient response, 
This corresponds to using a phase lead network in conventional design 
methods. Other poles and zeros near the origin tend to cancel each 
other and have little effect on the time response, Thus, a pole-zero 
pattern, such as shown in Figure 4,5 consisting of a single real zero 
and a pair of predominate poles should serve as the starting point in 
the synthesis of transient response, 

The quantity n, can be readily determined by considering the 
second term in Equation 4,10, Since the cosine term is always less 
than unity, the output will be certainly less than 5% of its final value 
value when the following equation is satisfied, 


n 
lz] Se 0,05 


spear Cs] Sica (4, 20) 
2 K = | 


PIS 


O 


Taking the logarithm of both sides yields 


67 








Figure 4,5. A basic pole-zero configuration for adjustment 
of transient response 
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Q(z.) 

eS - Log K “Piz ) 
Ho) oS) eee (4,21) 
Log |z, | 








A response with no output at zero time requires the coefficient 
of the second term in Equation 4,10 to equal the reciprocal of the 


cosine term with angle of 


[ ang Q(z.) - ang (z 0 - 1) - ang P'(z.) ] (4, 22) 


For the pole-zero pattern of the type shown in Figure 4.5, the angles 
are those shown in Figure 4,6 and generally the expression 4, 22 is 
about 180° and the coefficient term is nearly unity. In this case Equa- 


tion 4.21 can be replaced by an approximated but very useful result 


ae 1.3 (4, 23) 


log Iz. | 

The solution of Equation 4.23 generally yields a noninteger value for 
Nn. and the actual value is the upper integer value of the n, obtained, 
The locus of Zo for constant n, are circles on the z-plane with 
centers on the origin, The circles are shown in Figure 4,7 and can 
serve as starting point for locating the predominate poles, 


The number of oscillations may be calculated approximately as 


follows 
Nee . settling time (4, 24) 
period of damped oscillation 
or 
n 0. 
INS ee (4, 25) 
oo Mt 








B = -/QA,) + /Zo-! +/Plzo) 





Figure 4.6, Various important angles 
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Figure 4, 7. 


z 
oO 





circles of constant n 


a 


a2 


For control systems with steady-state error toa unit step input, 
the constant term in Equation 4,10 differs from unity by the extent 
of this error, The value of K may still be calculated since the amount 
of the error is given inthe specifications, However, in most practical 
cases, the error is negligibly small and the same procedure previously 
discussed may be followed. 

The error constants of interest are the position, velocity and 
acceleration constants, These were developed in section 2,4 in terms 
of G(z) and repeated below for convenience, 


K* = lim QG(z) 
P z—>l 


REY 8 ce ie eo lh aes 
Vv Al 
(ee | 
Mee Se ina ie Lee 
a 2 


il Z-—> | 


Introducing a change of variable 
eS Aen | (4, 26) 


The error coefficients become 


K* = lim G(@) (4, 27) 
P 80 
KS i ometey (4, 28) 
Vv af 
6 —+>0 
K* 2 1 ae G(@) (4, 29) 
= im 4,29 
= T° 6-70 


In terms of the variable © Equations 4,2 and 4.3 are 
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< 


Q(@) = (0 +1-4q,) (4, 30) 


k= 


— 


5 


P(@) = (o> lee) (4, 31) 


* 


i= 


~—_— 


From Equation 4,1 , the open looptransfer function of the sys- 


tem shown in Figure 4.1 with H(s) equal to unity is 


W 
K 1 (@ +1 - q,) 





G(6) = — (4, 32) 
fT (@+1-p,) -K 77 (@+1-4,) 
=| k=] 


Substitution of Equation 4,32 into Equation 4,27 and taking the limit 





yields 
w 
K Jf (1-4) 
kK * = _ (4, 33) 
P m Ww 
T (l-p)- K WT -q) 
i=l k=1 
Increasing the value of K will effect an increase of ae ees 


second approach would introduce a pole near the position 1,0 anda 
zero between it and the origin, 


Incorporating the substitution 


l 
eo (4, 35) 
k = Wk 


into Equations 4,30 and 4,31 yields 
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W 

Q(6) = 7T (6 + — (4, 36) 
i= k 
m 

Sekar Gs = ) (4, 37) 
1=] ; 


Introducing the change of variable and the substitution above in Equa- 


tion 4.18 gives 


i « 
. 


K = a (4, 38) 


8 


feato 
i 
— 


Thus it follows, for a system with unit feedback and no steady-state 


error toa position input, the closed-loop function becomes 


a. (b, @ + 1) 

Wie) ee 2 (4, 39) 
He (a, @ + 1) 
Sl 


The open-loop transfer function is then 


WwW 
TT (b, @ + 1) 
{se (4, 40) 
Ir WwW 
TT (a,8+1)- 77 (be +1) 
i=l k=] 


Expanding the above equation 
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WwW “f;pw ew w v/ 

ies bee [pbs] -2 bi] +++ + OTT Be 

Kz | 2g kz) k= k=! | 
G(@) eee ee 


Z SY; w m 2 ri 
m w em ee om BS ae sooo ae 
+b + 5 2.24 ice bq a t e, 7, i 


AY. ai 
c—/ Kal Rel 


ist 


Multiplying the above equation by @- and then setting 0 to zeroas 


indicated by Equation 4,28 yields the velocity error constant 





fans = 1 (4, 42) 
i m 
pays 
i=l = 
Ww m 
If yy b - 3 a is zero the velocity error constant is infinite 
k=] i=1 


T K* = ——— (4, 43) 





1/T 
K* = ee (4, 44) 
V m WwW 
\ 1 a l 
a0 > eam 
i=l Pi k=l] Ik 
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Z 
K* = $< (4,45) 
Z 2 


W 
a i ; yy ' 
l-q Be Ba 
k=1 k i= Pi 


The denominator of the above equation can be reduced (large K* ) 
by placing a zero near the point 1.0 . However, this will have the 
additional effect of increasing M and a . This will be illustrated 
by the following example, Consider a simple system with a configura- 


tion suchas shownin Figure 4.5. 


K (z - q)) 
W(z) = a (4, 46) 
(Z-ui tjv) (z - eee 


Placing these poles and zeros of the closed-loop function into 
Equation 4,46 , placing the expression to zero, and solving for qo 


which would yield an infinite velocity error constant gives 


Zz 
(1 = U4) + a 


ae 2 2 aa (4, 47) 
2(1 - u_) 
O 
Assigning values of uy = One. V2 = 0.4 results ina qo of 


0.5 , This location for the zero will yield an infinite velocity constant 
but the overshoot, M _, for this systemis 0,6. If the zero had 
been placed at the origin the value of M _ would be a more reasonable 
0.39 but the value of T K* would be unity. Thus it appears that 
compensating with one zero does not give independent control over 


peak overshoot and velocity error constant, 
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An arrangement which does permit adjustment of Ke without 
appreciable effect on the transient response makes use of a pole and 
a zero placed very close together in the form of a dipole and near to 
the point 1,0 as shown in Figure 4,8 . The zero at q; and the pole 
at p_ are placed sufficiently close together that the residue of the 
transfer function in the pole at Ps is made negligible, Although P,. 
and qd, are nearly equal, there can be a significant contribution to 


the denominator in Equation 4,44 because the inverse values 


and are large numbers | 
ed ae if 





Consideration of the Equation 4.45 shows that a dipole, but with 
the zero closer to the origin, may be used to raise the acceleration 
constant, Furthermore, additional zeros would be required for stabi - 
lization purposes giving rise to a typical pole-zero configuration shown 


in Figure 4.9 . 


4.3 The Required Number of Poles and Zeros 

The general location of the poles and zeros to achieve various 
requirements was shown in the previous section, However, at the start 
of the synthesis procedure it is necessary to know the minimum num- 
ber of poles and zeros that are needed by a given application, Although 
mathematically permissable, a practical compensator cannot possess 
poles at infinity. Consequently, the excess of poles over zeros for the 
compensated or desired open-loop function must be at least equal to 
the excess for the plant ae » oince the closed-loop function zeros 
are the same as those of the open-loop and the number of poles are the 
same for either open or closed loop function, it can be said that 


oe 








Figure 4,8, Pole-zero configuration for large velocity error 
constant 





Figure 4,9. Pole-zero configuration for large acceleration 
error constant 
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(number of poles - number of zeros) 
W(z) 
‘a (4, 48) 


(number of poles - number of zeros) Gz) 


It would be of considerable advantage to avoid the lengthy pro- 
cess of taking the z-transform from the Laplace transform in order 
to determine the number of zeros of the plant, In general this can be 
done by inspection as will be illustrated by considering the open-loop 
systems shown in Figure 4,10, When the continuous system as 
shown in Figure 4,10a has two poles in excess of the zeros, output 
will be available at the second sampling, l.e., n2=41, Since the value 
of the output can be determined at the sampling instant by dividing the 
denominator of G(z) into the numerator, an output at this sampling 
instant requires that the number of z-plane zero be only one less than 
the number of z-plane poles, The number of zeros equals the number 
of poles when the Laplace transform poles exceed the zeros by only 
one, Also, in these cases, a zero is found at the origin. If a hold 
network (1) is employed, as shownin Figure 4,10b, the same pro- 
cess of reasoning indicates that the number of poles exceeds that of 
zeros by one; however, a zero need not necessarily be at the origin, 
A plant including a delay, such as Figure 4,10c , of less than one 
sampling period will give rise to the same plane pole-zero differ- 
ence as discussed for Figure 4, 10a except the zero will not gener- 
ally be found at the origin, For every complete sampling period that 
the delay involves, poles are placed at the origin without any additional 
zeros, The overall transfer function for Figure 4, 10d is the product 


al (z) ne (z) and the difference between the poles and the zeros 
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(a) The plant is a continuous network 





(b) The plant includes a hold network 






F(s) / 


(c) The plant contains a pure delay 


C2) 
Els) “ E®@ man eS C(s) 


(d) Parts of the plant are separated by samplers 


Figure 4.10. Various types of plents 
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of each function is additive in the final transfer function, If the poles 
and zeros differ by one, then the over-all transfer function will have 
an excess of two poles over the zeros, This 1s in accordance with phy- 
sical expectation since the second network does not receive an input 


until the second sampling instant, 


4.4 Open- Loop Root Locus 

The method of this dissertation assumes the closed-loop function 
to be known leaving the open-loop function to be determined, Fora 
system suchas Figure 4.1 the open-loop function G(z) andthe 


closed-loop function W(z) are related by 


G(z) 


W(z) = 1 + G(z) 


(4, 49) 


Solving for G(z) gives 


W (z) 


G(z) = aha iE) (4, 50) 


Hence, the poles or roots of G(z) are the solution of the equation 
1 = W{z) = O (4. 51) 


Both functions are rational polynomials in z so that 


Ww 
K (z i qi.) 
- K Q(z) cz! 
W(z) = Py = =~ ari (4, 52) 
a (z - P;) 
i= 


r 

TT 

Koa (z - c,) 

k A(z) aan is 
B(z) m 

TU (Zz ve P,) 


j=) 


G(z) = (4, 53) 
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Equation 4.50 may be also writtenas 


k A(z) _ K Q(z) 

“Bi:) ~*~ Piy)- KQ&) (Aoe 
The gain constant are the same, the zeros of the open-loop function 
are identical with the zeros of the closed-loop function and the poles 


of the open-loop function are the value of z for which 


eh 


W(z) = (4, 55) 


The locus of roots consist of all points on the z-plane for which the 
plane of the closed-loopis 2hwm (h may be any integer or zero) 
that is, 


ang [ W(z) ] = O0+2hn (4, 56) 


In terms of Equation 4,52 the phase of W(z) is given by 


WwW Yn 
Q = y ang (z-4,) - » ang (z - p,) (4,57) 
k=] i: 


The precise position of the poles on the root locus requires that 


| Wiz) | = 1 (4, 58) 
or that 
w 
K fs |z a qi, | 
_ = |] (4, 59) 
[z e P, | 
Sl 


The first step in applying this graphical root finding method is 
to determine the locus of roots, An approximate root locus can 
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generally be found very quickly and simply by applying some simple 
rules concerning the geometry of the ‘Zero phase lines in the z-plane, 
After the approximate locus has been obtained the exact location of the 
roots can be determined by use of Equation 4.57 and 4,.59* . A typi- 
cal open-loop root locus is shown in Figure 4,11 

The rules and their proofs for the approximate construction of 
the open-loop root locus are now presented, 

Rule 1. ''The loci are continuous curves starting at a pole for 
zero gain and terminating at a zero at infinite gain, "' 

Rule 2. "The loci exist on a part of the real axis where an even 
number of poles plus zeros or none are found to the right," 

This may be easily verified by considerationof Equation 4,55 
and Figure 4,11. 

Rule 3, If the number of poles, w, exceeds the number of zeros, 
m , thereare w-m_ branches which goto infinity. Further, the 


directions of the infinite asymptotes with respect we the real axes are 
a 


2h 7, 


ae (4, 60) 
where n=0,1, 2, 3 

This rule is proven by reference to Figure 4,l¢a and b, The 
scale of the axes in Figure 4,12b is increased over that of Figure 4,l2a 
reducing the location of the poles and zeros of W(z) to almosta point, 
As point P moves further away from the origin the angles subtended 


by the zeros or poles are essentially the same, The sign of the angles 


from the zeros is apposite to the sign fromthe poles. Since the angles 


*For this operation a device called the Spirule has been developed by 
Spirule Company, 9728 El Venado, Whittier, Calif, 
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z-plane 
x closed-loop pole 


0 closed-loop zero 





Figure 4,11, Atypical open-loop root locus 
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(a) 





(b) 
Figure 4,12. Proof of rule three 
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from the excess poles are all the same, the angle of each pole is 


21 


a if the sumis equalto 2m, The following table will 


prove useful, 


TABLE 4,1 


Direction of Asymptotes 


w-m Directions 
l 0 
2 -180° 0 
3 = 1208 0 +120° 
4 -90° 0 +90° lee: 
Rule 4, "All asymptotes intersect ata point u, on the real 
axis," 
m Ww 
Golda 
i=l ° Re 
where u = (4, 61) 
- w-m 


The proof of this rule can be demonstrated by consideration of 


the transfer function of W(z) written in the form 


W(z) = m (4. 62) 


Dividing the numerator into the denominator to give 


K 
W(z) = (4. 63) 


,w-m n (b, : a,) ,7-m-1 t... 
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The locus is described by equating the denominator to -K, _ that is 


W-mM w-m-l 


Zz + (b,-a +. %2 = Terk (4, 64) 


ae 
By the theory of equations the sum of the roots of Equation 4, 64 is 
by - a 


asa polynomial of w-m_ degree with the sum of the roots given by 


; « When the magnitude of z is large Equation 4, 64 behaves 


es 71 


w-m™ is greater thantwo, This is no limitation as an inspection of 


being a constant which is independent of gain as long as 


Table 4,1 will show, With by - a, independent of gain, the sum of 


the root remains at one point for any value of gain, Hence, the asymp- 


totic lines radiate from the point on the real axis given by 


u = Jy as (4, 65) 
. w-m 
Rule 5, '!'The point u, at which the locus breaks away from 


the real axis is found by proper equating of the sum of reciprocal of 
the distances between poles and zeros and the point u, as explained 
by the following: "' 

The appreciation of this rule and its proof is demonstrated by 
Figure 4,13, Atrial point is assumed at xe and moved a small 


es 
vertical distance Av off the real axis, To be a point onthe root 


locus this point must satisfy the phase angle criterian, That is 


-9 + 06, + 0, +6, = 2hao (4, 66) 


All the angles are small and can be replaced by their tangents, hence 


87 








Figure 4,13, Proof of Rule Five 


[S 


Figure 4.14. Proof of Rule Six 


Ov, 


x1 


Os 
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— _ AY . fe 2 ioe S| 
u 


in ae + 180° = 360° (4, 67) 
- 1 ays faaiagy ea aa) Meee <6 


: + —— + _———— 2 ae ~~ (4, 68) 
Ven al ane e 4 e e z 


The solution of Equation 4,68 for the value of u. is generally 
accomplished by trial and error, For the case where complex poles 
or zeros exist, the expression needed to determine the break away 
point is more complicated. However, it has been shown (21) that unless 
the complex poles or zeros are close to the point of emergence their 
contribution is slight. 

Rule 6, "The angle at which a locus emanates from a complex 
pole or terminates upon a complex zero is found by summing the angles 
to all the other poles and zeros, When emanating from a pole the 
angles from the poles are taken negative, When terminating ona zero 
the angles from zeros are taken negative, "' 

Consider a test point P. shown in Figure 4,14 as onthe locus 
only a very small distance removed from the complex pole. Then the 
angle P. can be considered the angle of departure of the locus from 


the pole. Since Po is on the locus the sum of angles must apply. 
m Tc = 0 4.6 
Sy ae hy SO ee ese la): 6) (4. 69) 
solving for 


@, = [ 6, + 6, - (63 + ©, + 6) J (4. 70) 


As point P approaches the pole, the angles to the right of the 
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above equation become simply the angle subtended between the various 
poles and zeros and the complex pole in question, Obviously, when 


the angle of arrival at a zero is determined a change of sign is required, 


4.5 Closed-Loop Root Locus 

Application of the closed-loop root-locus assumes that the poles 
and zeros of the open-loop function are known and that the poles of the 
closed-loop function are to be determined, For a system such as Fig- 


ure 4,1 the open-loop and closed-loop function are related by 


W(z) = a (4, 71) 


Hence the poles or roots of W(z) are the solution of the equation 


As discussed in the previous Section, both functions are rational poly- 
nomials in z, Further, the zeros and gain factors of both functions 
are identical, Thus, the poles of the closed-loop function are the 


values of z for which 


Glo es 


h 


i 
= 
ie 
Wn 


(a 
The precise position of the poles or the root locus requires that 
| G(z) | = 1 (4, 74) 


As inthe case of the open-loop root locus, the first step in apply- 
ing this graphical root finding method is to determine the locus of roots, 
The approximate root locus are easily found by application of the rules 
which have been presented (6) (15) in the literature, <A typical closed- 
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loop root locus is shown in Figure 4,15, 


4.6 Additional Design Aids 


It has been tacitly assumed, before the techniques of this disser- 
tation are employed, the optimum system performance by simple gain 
adjustments was found short of the required specifications. 

This preliminary gain adjustment can be aided by Figures 4, 16 
and 4,17, These Figures relate the transient response quantities of 
M and a for second order systems with no position error to the 
location of the closed-loop complex poles, Figure 4,16 applies toa 
system with sampler only, while Figure 4,17 applies to a system incor- 
porating a zero order hold network (1) (12) . Only the first quadrant 
is shown since complex pole-zero patterns are symmetrical about the 
real axis and poles gn the second quadrant seldom give rise toa satis- 
factory response, The figures are entered on the circle correspond- 
ing to the product of s-plane open-loop pole and sampling period, These 
circles are actually the closed-loop root locus of the complex poles 
and proceeding counter clockwise corresponds to increased gain, On 
the circles are indicated the values of ,M_ and a that a second 
order system with a pole at that location will produce. For complex 
systems, the figures can be used on the basis of dominance to deter- 
mine if the various values of M_ and oo will satisfy the specifica- 
tions with sufficient margin to provide for the deleterious effect of the 
neglected poles and zeros, If this appears possible, a more exact 
study can be made using the closed-loop root locus, The figures may 
also be helpful for making the first selection of the predominate poles 
in synthesis process, 
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Figure 4,15. A typical closed-loop root locus 
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Figure 4,16, Values of M and a for second order systems 
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Figure 4.17, Values of M and n_ for second order system 
with hold P 
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The dipole required to raise the steady-state error constant may 
require such close spacing as to render praphical methods impractical, 
Applying the procedures of this thesis analytically would require the 
finding of roots of polynomials. Since the approximate values of the 
roots are already known, this is not as difficult a task as may first 
appear. However, an alternate method avoiding these difficulties works 
directly with the open-loop function, A pole anda zero are introduced 
into the open-loop function of proper relative values to give the desired 
increase of steady-state error constant as indicated by Equations 2.64, 
2.65 and 2,66, The expression for the output transform can be found 


as 
C(z) = a ee. ay (4, 75) 


where A(z) and B(z) are the numerator and denominator of the 
open-loop function while k is the gain, The actual output at the sam- 
pling instants are readily found by division, 

Unfortunately, this method leaves the designer in doubt as to 
whether the final output time response will be satisfactory. Actually 
the above procedure is not so involved that trial and error is out of 
order, Further, as an aid inthe initial selection it can be shown that 
the closer the loop pole of the dipole is to unity the smaller will be the 
residue of the output transfer function in the pole, 

This effect of open-loop pole location on transient response when 
employed to increase the steady-state velocity error coefficient will 
be illustrated as follows, For simplification let the original value of 


K* and the sampling period be both equal to unity. Let it be required 


v2 





that the value of Ke be raised by a factor of five. If the open-loop 
pole and zero are designated by a and q and the closed-loop pole by 


Pp then it follows directly from Equations 2.65 and 4.44 that 





ee Sn: | 
5 = (4, 76) 
and 
l 
eee (4, 77) 
[tap - 1 


Solving Equation 4.77 for p yields 


- 4 
(4, 78) 
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Bi er atria 
The introduction of a dipole into the closed-loop function affects 
the response in two ways. 
(1) The original terms are modified, but slightly since, asa 
result of the fraction (z-q)/(z-p) which is essentially unity at the 


Original poles, the residues at these poles are almost unchanged, 


12) Quays Gena poser st heraddicae 


p- p 
Consideration of this term will show that with movements of the 
dipole near the point 1.0 the coefficient of this term will change by 
a factor almost directly equalto (q-p)/(p-1) . Figure 4.18, 
employing the previous expressions, gives the value of this ratio as 


the open-loop pole is moved closer to unity in value, 
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Figure 4,18, The effect of the position of the open-loop pole 
on the output 
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CHAPTER V 


APPLICATION OF SYNTHESIS PROCEDURE 


5.1 Introduction 

Chapters III and IV are concerned chiefly with presenting the 
synthesis philosophy and the techniques necessary to implement the 
philosophy. Inthis chapter the actual process of synthesis is illus- 
trated by a variety of different specifications, Also discussed in 
detail is how the methods of this dissertation may be applied to sys- 
tems which do not have a sampler inthe error channel. The chapter 
is concluded by a discussion on synthesis of systems with multiple 


inputs, 


5.2 Complete Example of Synthesis Procedure 

This section presents an example of synthesis using the method 
of this dissertation, an example purposely chosen to illustrate the | 
basic elements of the procedure, Synthesis is initiated by determin- 
ation of the specifications which include the fixed part of the system 
and the time domain Benen ee characteristics, The specifications 
are as follows: 

(1) The plant has a sampler in the error channel operating with 


period of 0,1 second and a transfer function given by: 


K (s +19) 
s(s+3, 7)(s+8)(st+25) 


oan = (5.1) 


(2) The time domain specifications: 
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M < Pee is 

n < 4 
P 

n x 10 
Ss 

K = 00 
P 

K* ek 60 

Vv 


Approximating the plant as a second order system and employ- 
ie Figure 4,15 as explained in Section 4,6 shows that specifications 
cannot be satisfied by adjustment of gain and a compensator will be 
necessary. Consequently, the synthesis procedure outlined in Figure 
3.3 will be required, 

Two items which must be decided before selection of the closed- 
loop pole-zero configuration are the difference between the number of 
z-plane poles and zeros of the plant and how many of the poles and zeros 
of the plant are to be retained in the desired transfer function, For 
the first item it is seen that the plant has one sampler, no delays and 
a numerator three degrees less than the denominator, hence the z- 
transfer function of the plant will have one more pole than zero, For 
the second item, the retention of plant poles and zeros is largely tied 
up with consideration of the physical equipment, In this case an infin- 
ite value for So requires the pole at the origin to be retained and 
arbitrarily the pole at 3.7 will also be retained. 

With these constraints the simplest pole-zero configuration can 
be a pair of complex poles and one real zero, Usin? Figure 4.15 
again as a guide the first choice of pole-zero selection is made with 


poles at 0.5 + j0.45 andazeroat 0.3. Using the correlation 


theorems 
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n = 8 
S 


K* = 13 


Vv 


The transient response requirements are met with sufficient mar- 
gin to permit realization of a higher velocity error constant by intro- 
duction of a dipole, First, however, the open-loop poles must be ascer- 
tained, This is readily done by the open-loop root locus shown in Fig- 
ure 5,1, The gain is fixed by including one pole at point 1.0 and then 
the remaining pole is found at 0.64. Also shown in this figure are the 
construction lines used for measurements and angles to determine n 
and M., n. is estimated by reference to Figure 4.7, 

The s-plane pole of 3, 7 which is to be included in the final open- 
loop function corresponds to a z-plane pole of 0,691. Since the first 
try has produced a pole not too far removed from this location, a 
closed-loop root locus is now formed based on the zero at 0,3, and 
poles at 1 and 0,691 . This locus is shown in Figure 5,2 . Choos- 
ing on the locus poles at 0,55 +j0.46 gives the following time 


specifications 


With the gain determined, the open-loop function is 


591 (z - . 3) 
(z - l){z - . 691) 


G(z) = (52) 


In accordance with the discussion in Section 4,6 an open-loop 


pole is added at 0.99 andazeroat 0,95, The value of K* is 
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Figure 5,1, 


First selection for pole-zero configuration and 
open-loop root locus, 
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Figure 5,2, Second selection for pole-zero configuration and 
closed-loop locus, 
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then given by 


ay ae . 591 (1-595) ene een eee (5, 3) 
i (i. = . 99h coe: 


Further, realization of the open-loop function as a continuous function 
without time delays requires a zero at the origin. To retain the dif- 
ference of one between the number of poles and zeros an additional 
pole is added at 0.001. This pole was purposely placed very near 
the origin to cancel the effect that the zero may have on the transient 
response, The desired open-loop system is thus described by the 


transfer function 


ee) oe A(z) " po olrez (2 = 9a (a=. 3) (5, 4) 


B(z) (z- 1)(z-. 99)(z-. 691)(z-. 001) 


The inclusion of the additional poles and zeros will result ina 
shift of the complex poles and a change in the transient response, The 
new pole positions could be determined and the correlation theorems 
employed once again, Buta simpler method would be to form the out- 


put z-transform function C(z) toa unit step directly from G(z) as 





C(z) = Z x _ KAZ) (5.5) 
z- l K A(z) - B(z) 


By simple division the output at the sampling instant is in this case 


C(t) = .591 ui(t-. 1) + 1,087 ui (t-.2) + 1,33 u i (t-. 3) (5, 6) 
jee 34 w(t]. 4) 4 1.212 uo (tee Oo5 ww (t-. 6) 
O O O 
+ 0,966 ui(t-. 7) + 0,935 uo (t-. 8) + .955 uo(t-. 9) 
+ .995 u (t-1,0) 
O 
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The transient specifications are then seen to be 


M = 1,34 


Converting the z-transfer function to its equivalent Laplace trans- 
form with poles in the primary strip is the next step, This involves 
taking a partial-fraction expansion of 


G(z) 
Z, 
and multiplying through by z , converting term by term to Laplace 
transform with the aid of Table 2,2 , and, lastly, collecting the terms 


over a common denominator, These three steps take the form shown 


below, 
G(z) 6, 69 8z Bones ~9382z 7 . 2462 (5.7) 
i II z-.99 z- ,691 Be a ON! 
Grow 6,698 e Sele nee i , 246 (5, 8) 
Ss s+.l s+ 3.7 s - 6,9 
G(s) 20.95(s + .512)(s + 15,94) (5.9) 


s(s +t .1)(s + 3. 7)(s + 69) 


The compensation is now determined by dividing the derived 


transfer function, Equation 5,9, by the plant Equation 5,1, yielding 


G' (s) _ 1207 95(st. 512 lst’ )(stilos4)(stZ5) (5, 10) 


(s + .1)(s + 19)(s + 69) 
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The network is not realizable in its present form but the addition 
of a fourth pole far out on the negative axis, has a negligible effect on 


the transient response, and yields a realizable compensation function 


20. 95(st. 512)(st+8)(st+15, 94)(st25) 


2 (5, 11) 
(s t+ .1l)(s + 19)(s + 69)(s + 100) 


G(s) 
This function is realizable in terms of resistances and capacitors by 
standard network synthesis procedures (12) (23). 

The preceding paragraphs have attempted to illustrate the basic 
synthesis procedures, However, experience, personal preferences 
and considerations of particular systems may indicate deviations from 
the initial steps outlined above. As for example, if the entire plant 
must be retained as part of the final desired open-loop function it may 
be expedient to begin the design on that basis, This would correspond 
to entering the Flow chart (Figure 3, 3) initially from the plant posi- 
tion, This versatility of design procedure is the strong point of this 
synthesis method, This versatility is furnished by the two root locus 
techniques which relate the movements of poles in one function on the 
other and the correlation theorems which permit interpretations of the 


pole-zero configuration in terms of time response, 


5,3 Further Examples of System Design 

In this section three systems are considered, each with the same 
continuous fixed components involving one integration and an s-plane 
pole at 0,431, The first system includes only the required sampler, 
the second employs a sampler and holding network, while the third is 
characterized by a pure time delay of one half second, giving rise to 
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the following three transfer functions. 


G,_ (s) = —————— ; (5.02) 
Hl? s(s + . 431) 
ivecen> K 
G (sa Ko (5.13) 
ae S s(s + .431) 
K - .5s8 
G (s) = ——~———— x e °° (5,14) 
oe (sites oy 


In each case the sampling period is one second, 
Besides requiring that the above fixed elements be included in 
the final open-loop transfer function each system must also satisfy 


the time specifications, 


M < es 
n < 3 

Pp 
K = |) 

Vv 


With all systems required to meet the same time specifications 
at the sampling instants and having the same principle s-plane pole, 
the same closed-loop z-transform should be satisfactory for all three 
systems, A pole-zero configuration with one pair of poles at 0.518 


j0.55 andazeroat 0.15 meets the above specifications with 


>» 
Ag < 
iT 
—- WN 


Forming the output transform to a step function and dividing out 
gives the output time sequence as 
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the following three transfer functions, 


K 


G (s) = —_————— (512 
ne s(s + . 431) 
Ieewewe Ke 
G (s) = —— XK (5. 13) 
Be Ss s(s + .431) 
K - .5s 
Ce  —————— ee (5,14) 
ae s(s +. 431) 


In each case the sampling period is one second, 
Besides requiring that the above fixed elements be included in 
the final open-loop transfer function each system must also satisfy 


the time specifications, 


M < es 
n < 3 

Pp 
K = 1 

Vv 


With all systems required to meet the same time specifications 
at the sampling instants and having the same principle s-plane pole, 
the same closed-loop z-transform should be satisfactory for all three 
systems, A pole-zero configuration with one pair of poles at 05+ 


jO.55 andazero at 0,15 meets the above specifications w,ith 


M = oS 
n = 2. 56 
n = 3 

p 

K = I. 3 


Forming the output transform to a step function and dividing out 
gives the output time sequence as 
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alsin ee (Caps uo (t- 1) a ga ly 74) 72 ui (t- 2) tol 08 uy (t-3) 


to, 30 u (t-4) eval k Be (3; 3) 


a 


The reason for the difference between the predicted value of M_ and 
its actual value is due to the noninteger value for n as explained in 
Section 4.2. 

The closed-loop transfer function has a corresponding open- 


loop transfer function 


G(s) Sea ee (5,16) 


(Zee (zee 5) 


It should be noted that the plant pole has been accommodated 

since 
a+ Bol ace 

Although to this point the procedure has been the same for all 
three cases, the conversion to the Laplace transfer function, the next 
step in deriving the compensating network, is different for each case 
and will be illustrated in turn, 

(1) For the first system the plant is a continuous network whose 
z-transform is characterized by a zero at the origin, Adding this zero 
to Equation 5.15 anda pole 0,001 to maintain the proper difference 
between poles and zeros, the resulting transfer function becomes 


~ moo) z(ceen, 15) 
Gy gq (2) ST (5.17) 


(z - 1l)(z - .65)(z - . 001) 
Employing the methods described in the previous section yields 
the equivalent Laplace transfer function 


LON 





ee (en _ it o(s + 2, 82) (5. 18) 


s(s + .431)(s + 6.9) 
the transfer function of the compensator is then 


0,408(s + 2, 82) 
(s + 6.9) 


G(s) = (5,19) 


(2) Fora system including a zero orderhold network the trans- 


fer function of the open-loop system is the general form, 


x G. G. (z) (5720) 
Zz, 
The (Zee ee term corresponds to the transfer function of 
l-e ae , while the transfer function G_(s) of the pure integra- 


tion supplied by the hold network and the transfer function Go(s) of 
the continuous portion of the system result in the z-transform 


G. G*5 6 (s) inthe form of Equation 5,20 yields 


G5, (s) = a po (S220) 
z (7 = bh) (265) 
and then 
GHGetit| we 2 ee (5, 22) 
s 20 > 
(2 - 1) (z - . 65) 
The conversion of G. G59 (z) to its corresponding Laplace trans- 


form is similar to the previous example yielding 
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.449(s + 1.514) 


G (s)G (s) = (5, 23) 
‘ ae s (s + .431) 
hence 
or (s) = __8.98(s_ + 1,514) (5,24) 


s(s + .431)(s + 20) 


Where the pole at 20,0 has been added to enable physical reali- 


zation of the compensator 


0,0757(s + 1.514) 
(s + 20) 


G5. (s) (5, 25) 


(3) The third system includes a pure time delay of one-half the 
sampling period. Because of this delay a zero at the origin is not 
required inthe transfer function, Proceeding directly to the partial 


fraction expansion of Equation 5.16 gives 


G (2) Slee " ae 


—_ (5, 26) 
7, aoa Z => 


The conversion of each term on the right hand side of the above 
equation to the s-domain must be done through use of Table 2,3 , 


Modified z-transforms, with m equal to one-half, giving rise to 


CS) 
3d Z 1, 58 fi 1,15 (5. 27) 
Rae Ss s + , 43] 


Combining terms under a common denominator, adding a pole at twenty 


and making the necessary gain adjustment gives 
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- ,5s8 
8. 58(s + lee) S (5, 28) 


G (s) = 
ae sist JAG ene Or 


the required compensator is 


0.075(s + 1, 52) (5,29) 


G,_ (s) 
3¢ ie ow) 


A second series of three applications evolves when the value for 
Ke is increased to six while holding all other specifications the same, 


The open-loop z-transfer function satisfying these conditions is 


_ 65(z - 92) (z - 15) (5, 30) 


G g (z) = 
(z - l)(z - .98)(z - . 65) 


Employing previously discussed methods of this chapter the 


required open-loop systems are described by the transfer functions, 


1, 574(s +. 083)(s + 2, 88) (5, 31) 
Gig (s) Se ee 
s(s + .02)(s +. 431)(s + 6,9) 
Se 
ee (cee l-e ie es 9.477(s +. 8)(s + 1, 43) (5, 32) 
Ss s(s +.02)(s + .431)(s + 20) 
2 (G2 OLONG awa 2 ea 
6,6) 2 ee (5, 33) 


s(s t .02)(s + .431)(s + 20) 


The output time sequence to a unit step input for the above sys- 


tem is found to be 
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Cltijw= O65 ui (t-1) + 1,24 u (t-2) + 1. 48 uo (t-3) 


+ 1,37 u(t-4) + Bee (5, 34) 


Each of the six systems synthesized in this section were instru- 
mented on an analog computer. The computer outputs are given in 
Figures 5,3 through 5,8, It is to be noted that the output at sampling 
instant checks very closely with the calculated output. T#ble 5.1 com- 
pares the calculated and measured values for RY : 

Liberty was taken in Figures 5,5 and 5,8 where the curves are 
actually the input to the delay network displaced by half a sampling 
period, The imperfection of the computer time delay required this 
modification of the actual data, Justification for this procedure as 
well as detailed description of the computer study is given in Chapter 


Vi. 


5.4 Sampled Data System without a Sampler in the Error Channel 
The sampled data systems to be discussed in this section are 
characterized by a lack of sampling in the errorchannel, Figure 5,9 
is the block diagram of such a system, Reference to Table 2,2 gives 
as an expression for the output 
R G, (z) G, (z) 


C(z) = (5, 35) 
1+H G, G, (z) 


It is of particular interest that the transform of the input signal does 
not appear as a multiplicative factor in the numerator of the expres- 
sion for C(z). Consequently, it is not possible, as was done pre- 
viously, to divide out the transform of the input signal and forma 


closed-loop transfer function independent of the input, 
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Table 5,1, Comparison of calculated and measured 
velocity error constants 


System Calculated Measured 
sampler [258 ee [exe 
hold | Reaoy ta 1, 64 
delay 1,58 55 
sampler be 2 5.4 
hold Grae 2 at | 
delay oar 6-0 


sas 





Figure 5.9. 


H(s 


A sample-data control system 
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By the expediency of dividing both sides of Equation 5,35 bya 
z-transform of the input a closed-loop transfer function can be defined 


as 


RG {z7°G (z) / R(z) 
og) (5, 36) 


R 
1 + HG, G, (z) 





It should be noted that the transfer function applies only fora specific 
input and for this reason will be referred to asa limited closed-loop 
transfer function with designation Ww, (Zane 

In general, limited closed-loop transfer functions do not have the 
same zeros as the open-loop transfer functions even though the system 
provides a unity feedback path. Because of this characteristic it is 
not possible to employ the open-loop root locus for finding the poles of 
the open-loop functions nor are the expressions for the steady state 
error constants in terms of closed-loop poles and zeros applicable, 
However, the response to any particular input can be evaluated dir- 
ectly in the manner employed in Section 4,2 fora unit step and the 
response shaped as desired by proper selection of the closed-loop poles 
and zeros, The steady state error can be found by direct application 
of the final value theorem, Equation 2,55, For systems with zero 
steady state error toa step input, Equations 4,14 and 4.19 are 
directly applicable for synthesis to a prescribed transient response, 

With the limited closed-loop transfer function determined, the 
desired open-loop transfer function may be found by a method such as 
the following, Witha step input and setting H(s) to unity the remain- 


ing expressions defined in Equation 5,36 are polynomials in z_ with 


ly 


expressions 


RG) 2 (5. 37) 


a oral 2) 





R G, (z) = (5, 38) 
(z-1) b, (2) 
k, a, (z) 
et) = (5, 39) 
b, (z) 
Kea (z) 
G, G, (2) = —-—" (5, 40) 
by (z) by (z) 
= (z) _KQ (z) (5. 41) 
P (z) 


In Equations 5,38 and 5,40 , which are z-transforms of the product 
of two Laplace transforms, the convenient property that s-plane poles 
have a one to one image on the z-plane is utilized by factoring the 
denominator into the parts corresponding to each Laplace transform, 


Substitution of the above expression into Equation 5.36 yields 


la Kk ny 2,1 (2) k2a2(z) 


z) 


(5, 42) 
z ieee + kj), 22 (2) 


Generally, either G, (s) or G5 (s) are known in part or 
whole. Assuming that G, (s) is fixed, then the zeros of G, (s) 
are fixed by Equation 5,42 , after the closed-loop zgros and gain con- 


stant are selected, Consequently, the factors of b, (z) must be 
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chosen to satisfy Fquation 5,39 and also combine with the remaining 
open-loop poles and the open-loop zeros to form the poles of the closed- 
loop function with the exception of the pole at the origin introduced by 
the transform of the input, 

A first selection of the poles is aided by the consideration of the 
physical system and the type of compensation to be emploned: A phase 
lead network requires a pole near the origin, <A dipole for reduction 
of the steady state error requires a pole near the point 1,0. The 
total number of open-loop poles is one less than the number of closed- 
loop poles, The expression for A, (z) can be solved for by equating 
the numerators of each side of Equation 5,42 , A first trial expres- 
sion for G, (z) is then known and function G, (s) can be deter- 
mined by methods illustrated in the previous sections of this chapter, 
Then the expression for G, G, (z) may be formed and by root locus 
technique the closed-loop poles found. In general, these poles will 
not be those previously selected. However, the resulting poles may 
very well satisfy the specifications thereby completing the design, 

If this is not the case the process must be iterated. The precise 
method for determining Gy (s) may vary with different system 
arrangements, but with experience systems of this sort can be readily 


handled by application of the above principles, 


5.5 Systems with Multiple Inputs 

Because of the particular application or because of disturbances, 
all feedback control systems have multiple inputs. For example, Fig- 
ure 5,10 shows the block diagram of a sampled data servo system in 


which the output is corrupted bya spurious signal u(t) entering the 


Lis, 








Figure 5,10. Tandem single-loop system with secondary 
input 
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system at point A. Thus a feedback control system is characterized 
by n closed-loop transfer functions, the customary or primary trans- 
fer function relating output with the reference input and a number of 
secondary transfer functions for each secondary input to reference out- 
put. Whatever the source of the secondary input or wherever the 
point of entry to the system, satisfactory design involves the realiza- 
tion of suitable dynamic characteristics for each of the possible trans- 
missions, 

In Figure 5,10 G, (s) and G, (s) may represent transfer 
functions of fixed parts of the system while G. (s) is the transfer 
function of the compensator, Fora signal u(t) with transform U(z) 
the expression for the transform of the output is 


U G, (z) 
C(z) = ———————_—__—_——— (5.43) 


LG G, G, (z) 


All factors on the right side of the above Equation have been 
fixed by synthesis of the primary transfer function given in the speci- 
fications, No provision is available for modifying C(z) if the results 
of Equation 5,43 are unsatisfactory, More latitude is offered if the 
compensator transfer function G. (z) is split into tandem and shunt 
sections, This arrangement is illustrated by the system in Figure 
5,lla . Figure 5,11b shows the same system reduced toa single 
loop, The transfer functions of the blocks are related by the equations 


G(s) = G, (s) G, (s) (5, 44) 


12] 





U(s) 


+ t 
R{(s) SOs G, (s) = 1G, (s) srg cat 





G4 (s) 


wo 


(a) System with two input employing feedback and tandem compensator 


C(s) 


R(s) @ E(s) 4 


Gp (s) 





(b) Tandem single-loop system 


Figure 5,11. Two equivalent systems 
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G, (s) 
Gy (s) = ————— (5.45) 
1 + G> (s) G, (s) Gy, (s) 


The closed-loop z-transform transfer function for the system in Fig- 
ure 5,l1b is given by 
G. Gy) (z) 


(iz) = ————_-___— (5, 46) 
1 + G, Gy (z) 





R 


Figure 5,12a shows the system of the previous figure redrawn 
and the reference input omitted to emphasize the transfer character- 
istics of the system between secondary input and reference output, By 
block diagram manipulation the system can be reduced to that shown 
in Figure 5,12b where the transfer functions G. (s) and Gy (s) 
are the same as defined by Equations 5,44 and 5,45 ,. The transfer 
function of the output is piven by the expression 

UG (z) 


C(z) = —————_—_— (5,47) 
1 + G. Gy (z) 


Using the methods of the previous section the limited closed- 


loop transfer function for the secondary input is obtained as follows 


U G, (z) / U(z) 
<A (z) = a (5, 48) 
1 + G. Gy (z) 





Noting that the denominators of Equations 5,46 and 5,48 have a com- 


mon term permits as a more convenient expression for 


l2Z3 





U(s) 





(a) 


U(s) + C(s) 


G. te or. 


i 


(b) 


Figure 5,12, Two equivalent systems 
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i Ga ez) 
2 GC) = E: 2) | (5, 49) 
U(z) G, G, (z) 


Equation 5,49 places in evidence the possibility of satisfying 


independently specifications for both the primary and secondary trans- 








mission characteristics. However, U(z) and R(z) must be the 
same function if the equation is to be meaningful. Poles of © (z) 
and S (z) are identical with the exception of an additional pole 


at the origin of the expression for ae (z) . Consequently, the 
design freedom for meeting specification on as (z) is limited to 
choice of gain constant and zeros of iC; G, (z) . 

The choice of zeros and gain constant for meeting specifications 


are aided by the relationships given in Section 4,2 , With the design 








of the system for < (z) already completed and with the zeros 
of 5 (z) selected Equation 5.49 may be solved for the remain- 


ing unknown U G, (z) . Hence 








U G, ieft = | = (z) | = co | U(z) G. G, (z) (5,50) 


By procedures described in detail earlier in this chapter the transfer 
function for G,, (s) may be evaluated. The transfer function of the 


feedback compensator G, (s) is then given by 


G, (s) = eee! Bae 2 di | | <a (5,51) 
G, (s) G, (s) G, (s) 
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It is apparent that the poles of G, (s) and G, (s) , which are 
also poles of the open-loop function, are zeros of G, (s) . Conse- 
quently Gy (s) will in general require one or more poles at infinity. 
Hence, in practice G, (s) canonly be realized approximately by R- 
C circuits and high gain amplifiers. 

The expression for the transfer function G, (s) , given in 
terms of the previously determined desired open-loop transfer func- 
tion, Gy (s) , is 

G, (s) 


G, (s) = ————————- (5, 52) 


5.6 Suggestions for Future Work 

The field of sampled-data theory and application is still in its 
infancy. Many problems remainto be studied, One field that should 
prove very useful is the use of sampled-data theory in the analysis of 
continuous linear and nonlinear systems, In this case continuous func- 
tions would be represented by sampled approximations, The analysis 
of sampled-data systems with nonlinear components offers ample pos- 
sibilities for study. These areas of sampled-data theory are not dis- 
cussed in this dissertation, 

Of the problems which are a direct outgrowth of this dissertation 
three are of particular interest: first, extention of this synthesis pro- 
cedure for systems with non-periodic operation of the various sam- 
plers in the systems and consideration of wide sampling pulses which 


cannot be replaced by impulses, Second, further extention of the 
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procedure to control the shape of the response between the sampling 
instants, The third problem concerns the need to study the zeros of 
a z-transform to develope a method of predicting the location of zeros 


from the time function, 


Za 





CHAPTER VI 


DESCRIPTION OF THE ANALOG COMPUTER SET-UP 


An electronic analog computer study was undertaken to verify the 
synthesis methods of this dissertation, The computer was set up to 
correspond ona transfer function basis with the physical system, 

This is shown by Comparing Figure 6,1, where a block diagram ofa 
sample-data position servo is shown, and Figure 6,2 , where its com- 
puter analog is shown, The systemconsists of a sampler, phase lead 
network and a motor. In practice the latter device may be a commu- 
tator or any of a variety of analog to digital converters (23) (24). <A 
photograph of the entire computer layout is shown in Figure 6.3 . The 
important components of the set-up are identified on the figure. 

The computer set-up ona transfer function representation of a 
sampled data system served as more than just an aid to calculation, 
but as a direct analogy of the physical system and consequently a check 
on the analytical methods of this dissertation, Since the methods of 
this thesis are based on linear theory, the electronic analog computer 
simulation provides a more rigorous test thana physical Nyatcie No 
excuse is possible that nonlinearities obscured the results predicted 
by the theory. Also, the versatility of the computer permitted the 
testing of more possible system variations, 

As well as the system described above, a second system using a 
hold network and a third system incorporating a pure time delay as 
well as the sampler were tested, Each of the systems were designed 
to meet the same specifications and consequently the z-plane pole-zero 
configuration was the same, Two compensating networks were provided 
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Phase 
Lead 
Network 





= L6(S th.B2) 
Gis) $($ +-43/Xs +69) 


Figure 6,1. Simple sample-data servo system 
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for each system--the first included only a phase lead network to adjust 


for rise time and peak overshoot while the second added a dipole to 


Ag 


increase the velocity error constant. The details of the procedure for 
design of the compensators were explained in Section 5,3 . Also, 
shown in Figures 5.3 through 5.8 are copies of the output voltage ver- 
sus time curves secured from the computer set-up, These curves 
along with the tabulation of values for Kx in Table 5,1 represent 
the observed data of the computer study. The very close check 
between the analytical results and the observed results verified the 
methods of this dissertation and no further consideration of this mat- 
ter will be made inthis chapter, Rather, this chapter will limit its 
discussion to the new components which were developed to simulate 
the various system characteristics which are peculiar to a sample- 
data control system, The components referred to are (1) the Sam- 
pling unit, (2) the holding circuit, (3) the pure time delay, and (4) 
the input signal synchronizer circuit, 

(1) A front view of the Sampling unit (16) is shown on the left of 
the Figure 6,3 . Four sets of terminals are shown since actually this 
particular device furnished four sampling switches, A top view clearly 
exhibiting the various parts making up the device is available in Figure 
6.4 . From this figure it is seen thata D.C. motor* drives four ganged 
differential gears through a step down gear train, The output shaft of 
each of the four differential gears drives a split cam. Adjustment of 
the cam permits control over the switch closure time. Ford automo- 


bile distributer contacts are used for the actual switch, Each differ- 


*Electric Indicator Co., 1/25 HP, 115 volts, style ASC motor, 
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ential gear unit is furnished with a knob whieh permits displacing the 
angle alignment between its input and output shafts, Since the input 

of each differential is driven off the previous differential output (except 
for the first) this permits relative phasing between the closures of the 
four sampling switches, 

The Electronic Counter* shown in Figure 6,3 was used to mea- 
sure the sampling period and sampling duration, The sampling switch 
contacts were readily transferable from the computer circuit to the 
input of the counter in parallel with a battery and resistor circuit by 
means of the knife switch shown in front of the Sampling unit. This 
permitted control of the sampling period by adjustment of the potenti- 
ometer in series with the D.C, motor field. The electronic counter 
was adjusted to give time interval measurements to a tenth of a milli- 
second, Aone second sampling period was chosen and it was found 
that the motor speed was sufficiently steady to hold the sampling period 
variation to less than £5 milliseconds during a computer run, 

The closure time of the sampling switch was set for 59 milli- 
seconds, Employing Equation 2,2 indicates that the actual gain of 
the computer simulated system needed to be 16.95 times the calcu- 
lated gain in order that the pulse correspond in strength to an impulse, 
That the system would react to this pulse in approximately the same 
manner as to an impulse is assured since the principle time constant 
of the simulated system was 2,3 seconds, The sampling period of 
one second and the 59 milliseconds pulse were chosen as a compro- 


mise between the gain and impulse approximation requirements, 


*Hewlett Packard 522B Electronic Counter, 
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(2) The sampler plus the holding circuit was simulated by the 
circuit depicted in Figure 6.5 , It has the advantage of employing 
only one computer amplifier and the switch was furnished by the sam- 
pling unit, For proper operation the circuit functions as follows: 
when the switch closes the amplifier output voltage very quickly 
assumes the same value as the input (with reversed sign); when the 
switch opens the output voltage due to action of the capacitor will 
remain constant until the switch closes next time, It should be noted 
that the closure time of the switch should be ample to insure that the 
capacitor will charge up re its correct value of potential, With a 
switch closure time of 15 milliseconds suitable parameter values for 
R and CC were foundto be 0.01 megohms and 0.015 microfa- 
rods respectively. The preceeding values give a time constant of 0.15 
milliseconds, Since the sampling switch is closed for ten time con- 
stants it can be assumed that the difference between the input and out- 
put signal magnitudes is negligible. 

(3) Pure time delays or transportation delays can only be approx- 
imated by linear continuous networks, For this computer study the 
Padé fourth order approximant (12) was simulated by the computer 
circuit shown in Figure 6.6 (25) . The value of time delay [T was 
set for 0.5 seconds and the six operational amplifiers were secured 
from the separate Donner unit shown on the right of Figure 6.3 . In 
Figure 6,7 the output of this circuit to an input ramp signal is shown, 
It is significant that the output is not identically zero during the time 
delay but rather performs a low amplitude oscillation, The oscilla- 
tion is more pronounced if the input signal contains a large proportion 
of high frequencies, For this reason, the delay circuit was placed 
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Figure 6,5, 
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The sampler plus the holding circuit are 
simulated by this circuit 
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Figure 6,7. Input and output signals of the simulated 
time delay 
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at the output of the open-loop system permitting the system time con- 
stant to reduce the magnitude of the complementary frequency compon- 
ents introduced by the sampler. Nevertheless the output of the delay 
still contained an oscillation not present in the input signal. However, 
the system operates open circuit during this period and the signals do 
not interfere with the proper operation of the system since the delay 
assumes the proper output value by the time the sampling switch closes 
again, 

(4) Calculation of transient responses by analytical methods is 
usually performed with both the initial conditions and the input function 
zero at time zero, Fora sampled data system this requires input 
functions, such as ramps, to be first applied at a sampling instant, 
For a computer study this would require that the start of a computer 
run be synchronized with a sampling instant. The circuit depicted in 
Figure 6.8 was found to provide this synchronization, This circuit 
utilizes a second sampling switch A from the Sampling unit which must 
be synchronized with the switch used for sampling in the computer, 
The start and stop switch are both closed when a computer run is to 
be initiated. Then upon closure of the sampler contacts A , relay 
contact B closes a series circuit consisting of the coil and battery 
and thus locking the relay in the operate position while closure of relay 
contacts C starts onthe computer run. The Donner computer fur- 
nishes external operate contacts which were thus employed, The com- 


puter run is terminated by placing the stop switch in the open position, 
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Figure 6.8, Input signal synchronizer circuit 
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